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Preparation of aluminosilicate precursor by mechanochemical 
method from gibbsite-fumed silica mixtures 
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Abstract. Aluminosilicate (mullite) precursor was prepared by mechanochemical treatment of gibbsite and 
fumed silica mixtures. The effect of grinding on its structure and thermal behaviour was examined by ^^Al 
and ^‘"^Si MAS NMR, XRD, DTA-TG and FTIR. Thermal treatment of this precursor led to the ciystallization 
of mullite at about 1200°C via a spinel-phase. ^’AI and ^’'Si MAS NMR spectroscopies show lowering of 
mullitization temperature, which is associated with increased homogeneity of the precursor. Mechanochemical 
treatment of gibbsite and fumed silica mixtures resulted in the formation of a more homogeneous aluminosilicate 
precursor. 

Keywords. Mullite; mechanochemistry; NMR. 


1. Introduction 

The homogeneity of aluminosilicate (mullite) precursors 
depends largely on the processing method, i.e. how to 
mix, precipitate, hydrolyse or react Si 02 and AI 2 O 3 
components (Schneider et al 1994). As part of work to 
evaluate the mechanochemical effect for the preparation 
of aluminosilicate precursors, various alumina-silica mix¬ 
tures have been milled in a planetary pot mill (Temuujin 
et al 1997a-d). The homogeneity of these mixtures 
depends on water and hydroxyl content of the starting 
materials (Temuujin et al 1997b,c) and it was noted that 
the hydroxyls must be associated with the alumina starting 
material but not necessarily with the initial silica. 

Silica gel, fused silica and silicic acid were used as 
silica source, and found that the homogeneity of resulting 
precursors depended on the siliceous starting materials. 
In the present work, aluminosilicate (mullite) precursor 
was synthesized by mechanochemical activation of 
gibbsite and fumed silica mixtures and the product com¬ 
pared with those previously obtained under identical 
conditions using silica gel, fused silica and silicic acid 
(Temuujin et al 1997b,d). 

2. Experimental 

The starting materials used in the investigation were 
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gibbsite (Showa Denko Co. Ltd., Japan) and fumed silica 
(Aerosil 200, Aerosil Co. Ltd). The composition of 
the mixture was 75 ± 1 wt% alumina/25 ± 1 wt% silica, 
ground for 20 h in a planetary pot mill (La-PO.l, ITOH 
Co. Ltd) at room temperature. The rotation speed of 
planetary pot mill was 400 rpm and milling media and 
pot were made from silicon nitride. The weight ratio of 
milling media to powder was 20 : 1 . 

After grinding, samples were characterized by XRD 
(RIGAKU, Geigerflex with monochromatized CuKa 
radiation), FTIR (SHIMADZU, FTIR 8200 PC with the 
samples suspended in KBr discs) and DTA-TG 
(RIGAKU, Thermoplus TG 8120). ^""Si and ""Al MAS- 
NMR spectra were obtained at 11-7T using a Varian 
Unity 500 spectrometer and Doty probe spun at 10- 
12 kHz. The "®Si spectra were acquired using a 90° 
pulse of 6 |is and recycle time of 100 s, and were 
referenced to tetramethylsilane (TMS). The ""Al spectra 
were acquired using a 15° pulse of 1 |is and a recycle 
time of 1 s and were referenced to A 1 (H 20 )^. 

3. Results and discussion 

The DTA-TG measurements (not shown here), indicate 
that the intensity of the endothermic effect at 260°C 
due to dehydration of gibbsite decreases with grinding 
but does not disappear. A new exothermic peak which 
appears at about 960°C is related to the homogeneity 
of the resulting precursor, as previously described by 
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Figure 1. XRD pattern of ground mixtures and after calcining 
at different temperatures (■, gibbsite; •, spinel phase; unmarked 
peaks, mullite). 


Temuujin et al (1997a,b). Figure 1 shows XRD patterns 
of the mixtures ground for 20 h and calcined at different 
temperatures. In contrast with silica gel and silicic acid, 
the XRD intensities of the gibbsite reflections decrease 
with grinding and almost disappear. In the X-ray powder 
diffraction patterns of the samples, calcined for 2 h at 
1000°C, broad reflections of y-Al 203 or spinel-phase are 
observed. In our previous experiments we have shown 
that the DTA exothermic effect at about 980°C is related 
to the formation of spinel-phase into which is incorporated 
a small amount of silica from the aluminosilicate precursor 
(Temuujin et al 1997c). In the present precursor, mullite 
crystallizes at 1200°C without any other minor phases 
and its XRD intensity increases with subsequent heat 
treatments. The thermal behaviour of this precursor can 
be correlated with the structural changes introduced by 
grinding. 

Grinding for 20 h causes several OH stretching vibra¬ 
tions at 3500 cm“' to merge into a single broad band 
and OH bending vibrations around 1020 cm“' to disappear. 
Also the absorption band at 1100 cm”' due to the tetra¬ 
hedral Si-0 vibration shifts to 1020 cm"'. This shift to 
lower wave number can be explained by incorporation 
of Al^* into the Si'*^ sites (Temuujin et al 1997a,b; Okada 
and Otsuka 1986). The absorption band at 800 cm"', arising 
from three-dimensionally bonded SiO^ (Liao and Senna 
1993), also disappears on grinding. 

Figure 2 shows the ^’Al and ^‘^Si MAS NMR spectra 
of gibbsite, fumed silica and the ground mixture. 



Figure 2 . 11 -7 T and ~’A1 MAS-NMR spectra of starting materials and ground samples. 
Asterisks denote spinning side bands. 
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Unground gibbsite shows the octahedral resonance 
at 8-3 ppm, and unground fumed silica shows ^‘^Si reso¬ 
nance in the SiO, region (-111 ppm). In contrast with 
silica gel and silicic acid (Temuujin et al 1997a,d), 
fumed silica almost contains no hydrated Si-0 units. 
Grinding of mixture introduces new “’Al resonances at 
624 and about 32 4 ppm. The latter has been ascribed 
to pentacoordinated Al, or to distorted 4-coordinated Al 
associated with tricluster units which occur in mullite 
(Schmucker and Schneider 1996). The “‘^Si spectrum of 
ground mixture shows the appearance of a new strong 
peak at about - 83 ppm which reflects the presence of 
a significant proportion of the Si associated with alu¬ 
minosilicate units (Temuujin et al 1997c). The ^’Al and 
’‘’Si MAS NMR spectra of the ground mixture suggest 
that grinding for 20 h produces homogeneous structures 
in which the Si and Al are intimately combined, consistent 
with the crystallization behaviour of this precursor. 

These experiments suggest that fumed silica is more 
reactive to mechanochemical activation than silica gel, 
silicic acid and fused silica. Fumed silica, like fused 
silica, contains few hydroxyl groups, but the homogeneity 
of fumed silica + gibbsite mixtures as reflected by NMR 
measurements is much greater than fused silica + gibbsite 
mixtures (Temuujin et al 1997b). The present results 
suggest that the mechanochemical formation of amorphous 
aluminosilicate does not require the silica component to 
contain large concentration of hydroxyl groups. A more 
important factor may be the fineness and dispersion of 
the starting siliceous component (BET specific surface 
area of fumed silica is 200m’/g and fused silica is 
4-3 m’/g), which can beneficially influence the prepara¬ 
tion of aluminosilicate precursors by the mechanochemical 
method. 


4. Conclusions 

Mullite precursors were prepared by a mechanochemical 
method using gibbsite and fumed silica as starting 
materials. Judging by the homogeneity of the precursors, 
fumed silica is more reactive than silica gel, silicic acid 
and fused silica. These differences are attributed to the 
high degree of dispersion of the fumed silica (particle 
size less than several tens of nanometers), which is more 
completely adsorbed on the active sites of ground gibbsite, 
facilitating the formation of amorphous aluminosilicate. 
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Synthesis, characterization and application of a derivatized acid salt 
of tetravalent metal: o-chlorophenol anchored onto zirconium 
tungstate 
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Abstract. Tetra valent metal acid (TMA) salt zirconium tungstate (ZW) has been synthesized, folloveed by 
its deriv^tization using ortho chlorophenol (ZWoCP). ZWoCP has been characterized for elemental analysis, 
spectral analysis (FTIR) and thermal analysis (TGA). Its chemical stability has been assessed in various 
acidic, basic and organic media. Ion exchange capacity (I£C) has been determined and distribution behaviour 
towards several metal ions in different electrolyte concentrations studied and a few binary separations 
achieved. 

Keywords. Tetra valent metal acid salts; derivatized TMA salts. 


1. Introduction 

Tetravalent metal acid salts are obtained as amorphous 
and crystalline materials (Clearfield and Stynes 1964). 
These compounds have the general formula M(IV) 
(HX 04)2 • nH^O, where M(IV) = Zr, Ti, Th, Ce, Sn, etc 
and X = P, Mo, W, Sb, As, etc. The protons of the OH 
group can be exchanged for several metal ions and thus 
these materials are termed as inorganic ion exchangers. 

TMA salts can be derivatized by organic moieties 
bearing ionogenic groups such as -OH, -COOH, -SOjH, 
etc which also act as ion exchangers. These compounds 
are known as inorgano-organic ion exchangers or deri¬ 
vatized tetravalent metal acid (DTMA) salts (Alberti et 
al 1978, 1980, 1982; Kornyei 1982). 

Interest in DTMA salts stems from the fact that the 
inorganic backbone provides a stable and relatively inert 
anchor for the organic moiety. DTMA salts have the 
added advantage of both their counterparts in terms of 
thermal and chemical stability and ion exchange capacity. 
Derivatization of the inorganic ion exchangers by organic 
molecules depends on the nature of the inorganic matrix. 

In layered materials the organic molecules are incor¬ 
porated within the layers (Maya and Danis 1980; Costan- 
tino et al 1986), whereas in non layered materials the 
organic molecules get either anchored, sorbed on the 
surface or get encapsulated within the pores of the 
matrix. Preparation of organic derivatives of inorganic 
TMA salts has been studied (Alberti et al 1978, 1980, 
1982; Kornyei 1982) by intercalation of polar organic 
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molecules such as alcohols, glycols, acetone, urea and 
its derivatives or molecules of biopharmaceutical interest 
such as amino acids, heterocyclic bases, etc (Behrendt 
et al 1976; Costantino 1979; Maya and Danis 1980; 
Szirtes et al 1988) within the layers of crystalline a- 
zirconium phosphate (aZrP). 

Three oxygens of each phosphate group are bonded 
to three different zirconium atoms while the fourth 
oxygen bears the exchangeable proton that can be easily 
replaced by other cations. Compounds in which the =P-R 
or =P-OR groups (R = organic alkyl/aryl radicals such 
as H 2 O 3 P • CHjOH/HjOjP ■ C^H^) are reported (Alberti et 
al 1985: Alberti and Costantino 1991). 

Apart from the fundamental aspects, there is an 
increasing interest in the potential applications of inorgano- 
organic ion exchangers. According to the nature of R 
group, different applications in the field of chromato¬ 
graphic separation (Alberti et al 1989), catalysis (Alberti 
et al 1993a), photochemistry (Vermeulen and Thompson 
1992), protonic conduction (Casciola et al 1990; Alberti 
et al 1993b; Alberti and Casciola 1997) and radiochemical 
studies (Tomita Isao 1996), are possible. Organic 
derivatives are used as supports in gas chromatography, 
thin layer chromatography and molecular sieves (Alberti 
et al 1993a). The loading potential of the inorganic 
backbone is as high as two alkyl groups per zirconium 
atom, which is a distinct advantage over commercial 
reverse phase chromatography supports. Besides, inor¬ 
gano-organic materials are known to have a large se¬ 
lectivity for mono and multivalent cations (Costantino 
1979). 

Organic derivatives of layered aZrP can be termed as 



intercalation compounds or lamellar compounds. In case 
of yZrP however, derivatization takes place through a 
topotactic reaction which consists of the substitution of 
O3POH tetrahedra in a preformed layered phosphate by 
O3PR tetrahedra, whilst retaining the original layered 
structure (Alberti 1990; Alberti et al 1992, 1993b). Using 
this method, substitution of the OH by R groups (Alberti 
et al 1994) takes place as follows 

M(IV) (HPO,), + 2 RPO 3 H 3 ^ M(IV) (RP03)3 + 2 H 3 PO,. 

In the present endeavour, an attempt has been made 
to synthesize an inorgano-organic ion exchanger ZWoCP 
by treating an inorganic ion exchanger zirconium tungstate 
(ZW) reported by us earlier (Kale et al 1994) with an 
organic moiety ortho-ch]owpheno\ (oCP). 

ZWoCP has been characterized for elemental analysis, 
spectral analysis (FTTR) and thermal analysis (TGA). 
Ion exchange capacity (lEC) has been determined and 
the effect of heating on lEC studied. Distribution 
behaviour for several metal ions has been studied in 
various electrolyte concentrations, on the basis of which 
metal ion sepai'ations have been carried out. Binary 
separation of Co(II) from Cu(II), Zn(II) and Pb(II) has 
been achieved. 

2. Experimental 

2.1 Preparation of zirconium tungstate (ZW) 

ZW has been prepared by adding an aqueous solution 
of sodium tungstate to an aqueous solution of zirconium 
oxychloride till complete precipitation. The pH was 
adjusted to - 8 using very dilute ammonia followed by 
washing with conductivity water till the washings were 
chloride free. The precipitates were dried at ~ 40°C and 
sieved to the desired particle size (30-60 mesh). The 
sample was then converted to the hydrogen form by 
immersion in 1 M HCl, the acid being replaced inter¬ 
mittently. It was then washed repeatedly with conductivity 
water till free of chloride ion and finally dried at 40°C. 

2.2 Preparation of organic derivative (ZWoCP) 

ZWoCP was prepared by heating under reflux 1 g of the 
acid treated ZW with 20 ml of pure distilled o-chlo- 
rophenol. The excess of the organic reagent was decanted. 
This was followed by washing the sample with alcohol 
and finally drying it at ~ 40°C. The amount of oCP 
anchored onto ZW was estimated as the difference 
between the initial and final concentration by the bromi- 
nation method (Vogel 1989). 

2.3 Characterization 

Zirconium and tungsten in ZW and ZWoCP were esti¬ 
mated gravimetrically as zirconium oxide by cupferron 


method (Elving and Olson 1955) and as barium tungstate 
(Vogel 1978), respectively. Carbon and hydrogen in 
ZWoCP were determined using a Coleman analyser model 
33. Thermograms (TGA) were recorded on a Shimadzu 
DT-30 thermal analyser at a heating rate of 10°Cmin“'. 
FTIR spectra were recorded on a Perkin Elmer model 
1720X equipped with Epson Hi 80 printer/plotter. Sample 
preparation was done using KBr as reference with 1% 
sample concentration. 

Chemical stability of ZWoCP was checked in several 
mineral acids such as HCl, HNO 3 and H 2 SO 4 ; bases 
such as NaOH and KOH; and organic solvents such as 
alcohol, acetone and diethyl ether. 500 mg each of ZWoCP 
was taken in 50 ml of different solvents and kept for 
24 h and the changes in colour, nature and weight of 
the sample were noted. 

The ion exchange capacity (lEC) was determined by 
the column method (Nabi and Rao 1981) and the effect 
of heating on lEC was studied by heating the exchangers 
at various temperatures in the range 100°-500°C for 2 h 
each in a muffle furnace. 

2.4 Distribution studies 

Ion exchange behaviour for several metal ions such as 
transition metals (Zn^% Cu"'^, Ni^'*', Mn^"^, Co^"^); heavy 
metals (Pb^^ Bi^"", Hg^"^) and rare earths (La^^ Ce'^'^, 
Th'^'") was studied and their distribution coefficient {Kf) 
determined. is defined as 

(mlg"'). 

where, / is the initial amount of the metal ion in solution, 
F the final amount of the metal ion in solution, V the 
volume of the metal ion solution taken for analysis 
(ml) and W the weight of the exchanger taken for 
analysis (g). 

The ion exchanger (0-1 g) was equilibrated with a 
0001 M (20 ml) solution of the metal ion at room 
temperature for 24 h with intermittent shaking. The uptake 
of the metal ion was determined as a difference in the 
initial and the residual concentrations of the metal ion 
using EDTA as a titrant. Distribution studies were also 
carried out for all the above mentioned metal ions in 
different electrolyte media such as 0-001 M, 0-1 M, 0-4 M 
HNO 3 and 0 05 M, 0-5 M and 1 M NH 4 NO 3 solutions. 

2.5 Metal ion separations 

For binary separations, 5 g of the ion exchanger in 
form was taken in a glass column (30cmx 1-1 cm dia). 
The column was washed thoroughly with deionized water 
and the mixture to be separated was loaded on to it 
maintaining a flow rate of 0-5 ml min"'. In all the sepa- 
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rations, the cation concentration was maintained at 
O-OOIM concentration and 5 ml of each metal ion solution 
was taken. The separation of Co(II) from Cu(II), Zn(II) 
and Pb(II) was achieved by passing a suitable eluant 
and the metal ion in the effluant was determined quan¬ 
titatively by EDTA titrations. 

3. Results and discussion 

Chemical analysis indicates the composition of Zr : W 
to be 1 :2. Based on this and the TGA data, the proposed 
formula for ZW is ZrO, • 2 WO 4 • 9 H 2 O. The number of 
water molecules was determined using Alberti’s formula 
(Alberti and Torracca 1972). In case of ZWoCP, elemental 
analysis indicates, % C = 4-908; % H = 2-488 (instrumen¬ 
tal elemental analysis) and %C = 5-76; %H = 2-12 (bro- 
mination method). 

ZWoCP, though stable in HCl (11-3 N) and H^SO^ 
(36 N), was stable in HNO 3 only up to 0-5 N concen¬ 
tration. ZWoCP was found to be stable in bases such 
as NaOH and KOH up to 2 M concentration, as evidenced 
by no change in colour, form or weight of the samples 
used. ZWoCP was also found stable in ethanol, acetone 
and diethyl ether. 

FTIR spectra of ZW (figure la) shows broad bands 



in the ~3400cm"‘ region'characteristic of asymmetric 
and symmetric hydroxo-OH and aquo-OH stretching fre¬ 
quency. A sharp medium band at 1620 cm"' is attributed 
to (H-O-H) bending frequency. The FTIR spectrum of 
ZWoCP (figure lb) shows additional bands in the ~ 1580- 
1460 cm"' region characteristic of the C=C skeletal inplane 
vibrations. Bands in the ~ 1390-1330 cm"' region are 
due to the interactions of the -OH bending and C-H 
stretching frequencies. Weak bands in the ~ 1410- 
1310 cm"' region and at 1200 cm"' are due to C-O 
stretching frequency. 

The TGA of ZW indicates presence of hydrated water 
which is lost very slowly between 100°-180°C. TGA 
of ZWoCP indicates a weight loss at ~61°C due to 
traces of solvent used for washing the exchanger after 
refluxing. There is a gradual decrease in weight up to 
500°C indicative of the dissociation and decomposition 
of the organic moiety from the support and the conden¬ 
sation of the structural hydroxyl groups. 

The EEC of ZW at room temperature is 0-58 meq g"' 
whereas that of ZWoCP is 0-65 meq g"'. A higher value 
of lEC in case of ZWoCP is due to the additional 
sites of the ionogenic -OH groups in the benzene ring of 
the o-chlorophenol, besides the electron withdrawing phenyl 
ring also makes the phenolic (-OH) proton more labile. 



Figure 1. FTIR spectrum of a. ZW and b. ZWoCP. 


Table 1. Effect of heating on lEC of ZW and ZWoCP. 


Temperature 

(°C) 

Colour change 

Ion exchange capacity (meq - g ') 

ZW 

ZWoCP 

ZW 

ZWoCP 

RT 

White 

Yellowish brown 

0-58 

0-64 

100 

Light brown 

Brown 

0-40 

0-37 

200 

Light brown 

Black 

0-36 

0-21 

300 

White 

Black 

0-26 

0-48 

400 

White 

Dark brown 

0-24 

0-41 

500 

White 

White 

0-22 

0-30 
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Table 2. Effect of the electrolyte concentration on the values of metal ions for ZWoCP. 


Metal ion 

D.W. 

0-01 M 
HNO 3 

0-1 M 
HNO 3 

0-4 M 
HNO 3 

0-05 M 

NH 4 NO 3 

0-5 M 

NH 4 NO 3 

1 M 

NH 4 NO 3 

Cu^-" 

83-89 

6-78 

6-78 

NS 

83-90 

83-90 

83-90 

Zn--" 

67-20 

NS 

NS 

■ NS 

7-29 

NS 

NS 

Co^^ 

23-65 

NS 

NS 

NS 

10-22 

6-78 

5-04 

Mn^"^ 

78-92 

NS 

NS 

NS 

7-32 

2-04 

NS 

Ni^^ 

16-88 

NS 

NS 

NS 

15-37 

13-63 

NS 

Pb^-" 

99-21 

8-55 

6-78 

NS 

87-05 

90-82 

95-88 

Hg'" 

14-36 

11 -94 

6-78 

NS 

13-76 

17-84 

20-81 

Bi^"- 

1151-70 

240-43 

41-46 

NS 

64-00 

53-84 

27-59 

La^-” 

NS 

NS 

NS 

NS 

NS 

NS 

NS 

Ce'*^ 

NS 

NS 

NS 

NS 

NS 

NS 

NS 

Th^ 

67-20 

65-15 

60-23 

53-30 

45-41 

45-00 

40-00 


■•’'NS, No sorption 


Table 3. Metal ion separation for ZWoCP. 


Separation 

achieved 

Eluants 

Amount loaded 
(mg) 

Amount eluted 
(mg) 

% Efficiency 

Co(II)-Cu(Il) 

0-05 M NH 4 NO 3 (Co) 

0-2829 

0-2593 

91-66 


0-01 M HCl (Cu) 

0-3433 

0-1907 

55-00 

Co(II)-Zn(lI) 

0-01 M HNO 3 (Co) 

0-2829 

0-2574 

91-00 


0-05 M NH 4 NO 3 (Zn) 

1-6077 

1-3714 

85-30 

Co(II)-Pb(II) 

0-05 M NH 4 NO 3 (Co) 

0-2829 

0-2652 

93-00 


0-01 M HNO 3 (Pb) 

1-0774 

1-0567 

98-00 


The results of the effect of heating on lEC are sum¬ 
marized in table 1. For ZW, the lEC decreases on 
heating. This is due to the condensation of the structural 
hydroxyl groups bearing the exchangeable protons, at 
higher temperatures. Hovt^ever in case of ZWoCP, the 
lEC decreases up to 200°C and then increases for the 
sample heated at 300°C. Again for the samples heated 
at 400° and 500°C, the lEC decreases. This may be due 
to the fact that on heating in the range 200°-300°C, the 
organic moiety decomposes leading to the formation of 
active carbon as evidenced by the change in colour of 
the sample to black. This is responsible for an increase 
in the lEC and is later lost as COj as the heating 
temperature increases. This fact is also evident from the 
IR spectra of heated samples of ZW and ZWoCP. 

The distribution coefficient {K^ calculated for the 
metal ions in different electrolyte concentrations for 
ZWoCP are summarized in table 2. On the basis of the 
selective affinity for particular metal ions in the presence 
of electrolytes, the metal ion Co(II) has been separated 
from Cu(II), Zn(II). and Pb(II). Co(II) has been eluted 
using 0-05 M NH 4 NO 3 in all the separations except for 
Co(II)-Zn(n) separation,,where it is eluted using 0-01 M 
HNO3. Cu(II), Zn(n) and Pb(II) were eluted using 0-01 
M HCl, 0-01 and 0-4 M HNO 3 solutions, respectively. 
The results of the separation study are summarized in 
table 3. 


4. Conclusion 

ZWoCP exhibits good ion exchange capacity and is 
stable up to a fairly high temperature, besides exhibiting 
stability in different acidic, basic and organic media. Its 
selectivity for particular metal ions enables it to act as 
an efficient ion exchanger for the separation of Co(II) 
from Cu(II), Zn(II) and Pb(II) with -93% efficiency. 
The DTMA ZWoCP thus exhibits the characteristics of 
a promising ion exchanger. 
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Abstract. A simple extension of Beale’s and Lehmann’s models for the formation of porous silicon layer 
on p-type silicon is proposed with a view to explain the experimental conditions necessary for obtaining 
either uniform vertical pores or non-uniform pore branching, as desired. A uniformity parameter is defined 
and correlated with the measured porosity. The dependence of the porosity and the uniformity factor with 
the various formation parameters of porous layer are studied experimentally and explained qualitatively. 

Keywords. Porous silicon; p-type Si. 


1. Introduction 

The discovery of room-temperature photoluminescence 
(PL) of porous silicon has initiated an extensive research 
on both theoretical and experimental aspects especially 
because of its strong potential in its use in optoelectronic 
and photonic applications. 

Porous silicon (PS) is formed by anodic etching of 
crystalline silicon in hydrofluoric acid. The specific struc¬ 
ture of PS depends on various parameters like doping 
type and level of c-Si, HF concentration, formation 
current density, presence of illumination, surface mor¬ 
phology etc (Smith and Collins 1992). Depending upon 
the various parameters during formation, PS may be 
looked upon as either quantum wires or quantum sponges 
(Koshida et al 1992; Marusak et al 1993; Feng and Tsu 
1994b; Jain et al 1994; Kocka et al 1996). 

For optoelectronic applications, however, the near¬ 
vertical growth of pores and consequently that of the 
silicon pillars should be a pre-requisite to ensure repro¬ 
ducibility as well as isolation of the devices formed on 
a single wafer. 

In this paper, we report a detailed study on p-iype 
silicon to examine the growth of PS under various 
experimental conditions with a view to find out the 
favourable formation conditions for near-vertical growth 
of silicon pillars. The formation model proposed by 
Lehmann is extended in this paper introducing the concept 
of a uniformity parameter (Lehmann 1993). 

2. Modelling of vertical PS growth 

The basic prerequisites for electrochemical pore formation 
*Author for correspondence 


are passive state of pore walls and active state of pore 
tips (Feng and Tsu 1994a). Dissolution of silicon in HF 
occurs if holes are present in the surface. Thus, a surface 
depleted of holes is passivated. Hole depletion occurs 
if any hole reaching the solid-electrolyte interface is 
immediately consumed in the dissolution reaction. This 
is the case when the chemical reaction is not limited 
by the mass transfer in the electrolyte—a condition 
satisfied when the formation current density is below a 
critical current density value 

However, when the reaction is limited by mass 

transfer of ions in the electrolyte. The dissolution of Si 
is controlled by the supply of reactants through the 
diffusion layer in the electrolyte, leading to a build up 
of surface charge in silicon due to accumulation of holes 
and resulting into electropolishing (Turner 1958). 

In equilibrium macropore formation, a steady-state 
condition is assumed at the pore tips, i.e. the charge 
transfer due to holes exactly balances the availability of 
ions in the electrolyte (John and Singh 1995). In such 
a case, all holes will be consumed at the pore tips itself 
so that the pores grow vertically downwards with a 
constant cross-section. This is characterized by critical 
current density 

Thus, if the localized current density, at any place 
is greater than then the steady-state condition for 
lateral electropolishing will be achieved when the cross- 
sectional area to increases to a value Ap^ so that 
= This concept is realized at the initial stage of 
formation of the porous structure. The initiation of PS 
structure probably occurs at the locations of random 
inhomogeneities in the path of current flow through the 
interface (Beale et al 1985). 

In regions of the inhomogeneities where the local 
current density is initially high so that J^.^ > , localized 
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lateral electropolishing should occur until becomes 
equal to 

If the total number of inhomogeneities (nucleation 
centres) at the surface is then the current flow 
per m^) occurs through an effective area of so 

that 

iNAJJ=JA, (1) 

A is the total cross-sectional area of the same and J the 
measured current density. 

If the growth of pores were uniform, then the porous 
silicon thickness would be given by 


U'l' 


of This happens because of the fact that increases 
in proportion of J for a constant 

If, however, the number of nucleation centres (A^p,) is 
very large, then the maximum limit to pore wall thickness 
is set by the pinch-off effect (Lehmann 1993) due to 
overlap of adjacent depletion layers so that the porous 
column dimension is twice the depletion width (AVK) 
(figure 1). 

Under these circumstances, we have 


P..= 


and 




J 


( 6 ) 




ne 


( 2 ) 




A 

(d + 2A Wf' 


where n is the number of atoms required for dissolution 
of each silicon atom and p the Si density and n, reported 
to be of the order of 2 (Beale et al 1985). 

The porosity of vertically grown PS (no branching) 
will therefore be given by 


p Kp _ J 

" Ah^p A J^' 


It can be seen from (l)-(3) that 

“ nep 
and 

Ph = - 

“ “ nep 


(4) 

(5) 


Thus, when PS growth is vertical, the porosity is a 
function of J and (in turn, dependent on HF con¬ 
centration). However, the thickness of PS layer is inde¬ 
pendent of the current density and is only a function 



since the number of pores per unit area is \l{d + 2AWf. 

Thus, in this case, the porosity becomes a function 
of the pore diameter and the sample resistivity in the 
form given by 


d- 


lAW 


iP 


J 


(7) 


and the formation of pores will be possible only if 
A . .1 


N <• 

4AW' 


i.e. < 


4AW^ 


( 8 ) 


It is quite obvious that heavily doped p-type silicon 
will give rise to a large number of nucleation centres 
in the vicinity of doping impurities (Beale et al 1985), 
since current flow can occur in such samples by tunneling 
and a larger fraction of area is then always amenable 
to offer the path of current flow. Due to the depletion 
layer width overlap and pinch-off effect, the columnar 
growth of the pores is highly likely in the case (Beale 
et al 1985; Lehmann 1993). 

Equation (7) is indeed obeyed in the case of heavily 
doped layers as can be seen by inserting the results of 
Beale et al (1985) into (7). The number of nucleation 
centres necessary for vertical pores is given by 


N. 


- 


1/2 

2 

rjf_\ 


-2AW 

= 



N 


J 

- 

1 ’’•V 

- 


1 ^ J 


(9) 


Thus, the number of nucleation centres per unit area 
(”psc) vertical PS growth (from (9) and (8)) is given 
by 


Figure 1. Schematic diagram of uniform pores where d is 
pore diameter and 2AW is the separation between two con¬ 
secutive pores. 



2AW 


( 10 ) 
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and the pore dimension d is given by 
2{AW) 

Thus the pore dimension is found to be strongly dependent 
on the doping concentration. 

It is obvious that if is less than then the pore 
walls will not be passivated by the pinch-off effect and 
the branching of the pores will be a possible phenomenon. 

The critical nucleation centre density decreases 
and thus approaches towards the actual centre density 
as J approaches and ‘AW decreases (i.e for 
lower resistivity) as tong as the actual is not dependent 
on either J or AW. 

3. Pore branching and uniformity 

The morphology of pores formed depends on a number 
of parameters such as the surface roughness, doping type 
and density, anodizing current and voltage, pH and 
concentration of HF solution etc (Feng and Tsu 1994c). 
Depending upon these parameters, the pores may have 
fir-tree-shaped structures with interconnected branches 
underneath or well-defined unbranched vertical wells. In 
general, unless the surface is very polished and inten¬ 
tionally patterned, doping density is high and the ano¬ 
dization potential, HF concentration are of right 
magnitude, the pores exhibit branching and non-uniform¬ 
ity (Lehmann 1993). 

The measured porosity (P,^) of PS will be different 
from that for uniform growth due to pore branching and 
can be expressed as 

p mass of silicon dissolved 
mass of the porous layer 


fJAA 


ne 

\ J 


X m. 


m f 


Jtni. 


nep 


h 


( 11 ) 


where is the thickness of the porous layer in the 
direction of external current density. Combining (5) and 
( 11 ), we have 


Ll, 


where u is defined to be a measure of uniformity of 
the formed porous layer. In the case of a vertical pore 
growth, we have u = \. The uniformity factor 

(u) can be estimated from the relation 




- . (L 

P P 

* m m 

The measured porosity (P,„) for a given J will be highe 


due to pore branching than that for a uniform PS growth. 
At lower current density holes are depleted 

from the silicon surface resulting into a preferential 
transport of holes from the bulk to the surface irregu¬ 
larities. This implies that higher current density will lead 
to lower non-uniformity (John and Singh 1995). 

4. Experimental 

Silicon wafers of / 7 -type and (100) orientation and 7 cm^ 
area, having resistivity of 1-2 Q-cm were anodized in 
a cell specially developed for the purpose. The wafers 
were texturized by standard anisotropic etching process 
that leads to the formation of pyramidal structures with 
exposed ( 111 ) planes. 

The wafer actually acts as a seal between the front 
and rear regions. The front region of the cell was filled 
with HF: CH3OH, while the rear portion was immersed 
in KCl solution. The back contact metallization was done 
by screen printing of silver aluminium paste and its 
subsequent firing. The textured samples were anodized 
by a constant current supply (Bandopadhyay et al 1996). 

The gravimetric method was employed to determine 
the porosity of the samples (Amato 1995). Samples, 
before and after anodization, were weighed by the micro 
balance model No ADN200W. Then, the porous layers 
were removed by dipping the samples in NaOH solution 
and the wafers were again weighed. 

The porosity was then determined by using the relation 

m, — m. 

p _ * _£ 

where m, is the mass of the sample before anodization, 
the mass of the sample after anodization and W 3 the 
mass of the sample after complete dissolution of the 
porous layer. 

SEM techniques have been used to provide some of 
the detailed direct information about the thickness of 
the porous layers and the branching of pores in these 
layers. 

5. Results and discussion 

5.1 Role of formation current density 

Figure 2 shows the variation of measured porosity (P^) 
as a function of current density (7) for a fixed HF 
concentration (mixed in 1 : 1 proportion with CH3OH) 
and fixed time of anodization. The porosity (P^) increases 
with increasing J. 

The thickness (/ij of the porous layer has been 
estimated from the mass of the porous layer and is given 
by 


h 


m 


mj - m3 

Ap 
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Figure 3 and table 1 show the variation of estimated 
with the current density (7). It is observed that 
increases with increase in J. Figures 4 and 5 show two 
typical scanning electron micrographs of porous silicon 
layers at 45° tilt angle. It is seen from these figures 
and table 1 that the thickness of porous layers as 
measured by SEM technique agree fairly well with those 
measured by gravimetric methods. 

It has been seen from the simple model presented 



Figure 2. Variation of measured porosity with cieliine 

current density (7) when other etching parameters are constant. 



Figure 3. Variation of porous layer thickness with etching 
current density. 


here that for a given passed charge {JAt) more branching 
will lead to a lower PS layer thickness and consequently 
higher porosity than that would be observed for uniform 
vertical pore growth. Thus, it appears from this study 
that as the fractional increase of P (AP /P ) is much 
less than the fractional increase of J (AJ/J) and since 
increases with increasing J, the average thickness 
(/^n,) approaches towards the uniform growth thickness 
(hj with increasing J. 

The degree of uniformity is examined by investigating 
the variation of the uniformity (u) with the formation 
current density (J). Figure 6 shows that with increasing 
J, u is found to increase almost linearly with J. Table 
1 shows the different values of P^, u and as 
functions of current density. It is seen from table 1 that 
ii lies in the range 0-04-0-17 for the current densities 
in the range of 10-70mA/cm^ This can be explained 
from the consideration that as the critical density of 



Figure 4. SEM of porous silicon formed under anodic etching 
by J=10i'nA/cm“ in 48% HF for 5 min, estimated 
=4-25 



Figure 5. SEM of porous silicon formed under anodic etching 
by 7 = 50 mA/cm* in 48% HF for 5 min, estimated li^ = 9-6 jum. 
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Table 1. Values of measured porosity, PS layer thickness for non uniform (h^) 
and 100% uniform {hj and uniformity (m) for different etching current densities. 


J 

(mA/cm") 

(%) 

(pm) 

(from gravimetric 
method) 

K 

m 

(pm) 

(by SEM tech.) 

u 

(%) 

K 

(pm) 

10 

48-7 

4-7 

4-25 

3-64 

129 

30 

78-2 

7-6 

7-8 

6-8 

112 

40 

83-1 

8-3 

8-6 

8-5 

98 

50 

83-3 

9-9 

9-6 

10-6 

93 

60 

86-8 

14-0 

- 

12-3 

114 

70 

89-4 

151 

- 

13-9 

108 



Figure 6. Variation of uniformity (u%) with etching current 
density. 

nucleation centres decreases with J almost linearly 
according to (10), an increase in J will favour the vertical 
pore growth since the actual density of nucleation centres 
approaches with increasing current density. 

It is apparent from the present study that in the range 
of current densities considered here (10-70 mA/cm^), the 
uniformity factor is not very high (<017). The uniformity 
factor is increased by increasing J, but it is likely that at 
high current density the pore tips lead to increased oxide 
formation (Feng and Tsu 1994c). Thus, it seems that it is 
preferable to either increase by pattern formation on the 
surface and creating the required value of ~ n^^ and/or 
lowering the value of by increasing AW or in other 
words, decreasing the doping concentration of the sample. 

The fact that the anodizing current mainly flows through 
the pores and not through the whole area (A) of the 
sample was examined by calculating the number of 
electrons required for the removal of every silicon atom 
estimated from the weight change of the wafer during 

annHi'/iirinn nc ■fnllnu/c 


Table 2. Number of electrons (n) required for dissolution of 
one silicon atom during electrochemical etching of p-type Si 
for different values of etching current density (J). 


J 

(mA/cm^) 

m, 

(g) 

/«2 

(g) 

(g) n 

10 

1 0251 

1-0214 

1-7 

30 

0-9987 

0-989 

1-9 

40 

1-1263 

1-115 

4-656x 10-23 2-2 

50 

1-0431 

1-0296 

2-3 

60 

0-9769 

0-9571 

1-9 

70 

1-081 

1-059 

1-9 


= — (13) 

where mj = mass of each Si atom. 

It is found that n is of the order of 2 indicating that 
the flow is essentially through the pores and on an 
average, two electrons are required for each silicon atom 
removal. It was, however, observed (table 2) that n 
increases in general with increasing J and values of n 
were found to lie in the range 1 •7-2-3. 

5.2 Role of methanol content in HF electrolyte 

The role of methanol in the formation of PS layer has 
been examined by anodization of silicon in the same 
HF concentration with different proportions of methanol. 

It is seen from table 3 and figures 7 and 8 that with 
increasing methanol proportions in HF, the average PS 
layer thickness and the uniformity factor (u) decreases 
for a given current density. This is possibly because of 
the fact that methanol acting as a wetting agent removes 
the hydrogen bubbles from the pore walls leading to the 
branching of pores. 

The pore formation rate (average thickness per unit 
time) is found to follow the relation 

(14) 

as suggested by Lehmann et al (Feng and Tsu 1994a) 
for 1:1 HF (48%)-methanol electrolyte. 

5.3 Role of HF concentration 

Fiaiirp.Q Q nnd 10 chnw thp nf norr»<iitv nnH 



found that porosity increases linearly with HF concen¬ 
tration, whereas the uniformity factor decreases. 

This can be explained with reference to the change 
in As HF concentration increases, the value of 
also increases leading to the increase of the critical 
nucleation centre density needed for near vertical growth. 
This leads to lower values of with increasing 


observation is further corroborated by reflectance studies 
of PS (to be communicated separately). 

The fact that increase of HF concentration leads to 
less uniformity and consequently more pore branching 
can be observed directly by examining the SEM mi¬ 
crographs taken from a direction perpendicular to the 
thickness of porous layer (figures 11 and 12). The black 


Table 3. Values of measured porosity (P^), PS layer thickness 
(h^) and uniformity (w) for different methanol proportions in 
the electrolyte. 





u 

CH 3 OH: HF 

(%) 

(ttm) 

(%) 

1 :1 

78-0 

6-3 

6-8 

I :2 

75-0 

74 

7-1 

1 :3 

744 

8-2 

7-2 

1 :5 

724 

8-9 

1-2 

0:1 

51-5 

9-9 

10-3 

Table 4. Values of measured porosity (Pj^), PS layer thickness 
(h^) and uniformity (u) for different HF concentrations. 

C 



u 

(wt%) 

(%) 

(wm) 

(%) 

48 

83 

2-5 

2-1 

36 

74 

2-1 

3-2 

30 

77 

1-3 

3-7 

24 

63 

2-3 

5-6 

18 

61 

20 

7-6 



0 12 3 4 5 6 

HF methanol 



HP ; methanol 


Figure 8 . Variation of uniformity («%) with HF, methanol 
ratio in the electrolyte. 



HF concentration 


Figure 7. Variation of PS layer thickness (h^) with HF, Figure 9. Variation of measured porosity iP^%) with HF 
methanol ratio in the electrolyte. concentration (wt%). 
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HF concentration 

Figure 10. Variation of uniformity (u%) with HF concentration 
(wt%). 



Figure 11. SEM micrograph of porous silicon (HF = 24%, 
7=10mA/cm^) showing branching of pores. 

holes in these micrographs show pores formed in the 
direction perpendicular to the thickness of porous layers 
and thus these indicate the branching of pores. Figures 
11 and 12 clearly show that more branching occurs at 
higher HF concentration. 

6. Conclusions 

A simple model based on Beale’s and Lehmann’s ideas, 
for uniform porous layer formation on p-S\ has been 
developed correlating the various growth parameters of 
vertical and uniform pores with the pore formation 
parameters. A concept of non-uniformity has been 



Figure 12. SEM micrograph of porous silicon (HF = 48%, 
y=10mA/cm^) showing branching of pores. 

introduced to characterize the quality of PS layer in 
respect of pore-branching and near-vertical growth of 
pores. A factor called ‘uniformity factor’ has been defined 
and its dependence on experimental and physical 
parameters of porous layer studied. The simple model 
developed here has been utilized to explain the different 
experimental observations. 

With increase in J and/or decrease in HF concentration, 
the uniformity factor increases leading to better vertical 
growth. With more methanol proportion in HF electrolyte, 
the uniformity factor decreases giving more branching. 
The results can be explained with reference to the critical 
density of nucleation centres on the silicon surface. 
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in a monotectic system 
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Abstract. The phase diagram of a binary organic system involving diphenyl and succinonitrile shows the 
formation of a eutectic (0-968 mole fraction of succinonitrile) and a monotectic (0-074 mole fraction of 
succinonitrile) with a large miscibility gap in the system, the upper consolute temperature being 53-5°C 
above the monotectic horizontal. From the enthalpy of fusion of the pure components, the eutectic and the 
monotectic, determined by the DSC method, the enthalpy of mixing, Jackson’s roughness parameter, interfacial 
energy, size of the critical nucleus and excess thermodynamic functions were calculated. 

Keywords. Phase diagram; monotectic system; thermochemistry; organic eutectic; organic monotectic. 


1. Introduction 

The fundamental understanding of solidification (Elliott 
1983; Herlach et al 1993; Glazer 1995; Rzyman et al 
1996) and properties of polyphase alloys has become a 
subject of potential investigation during the past two 
decades. Due to low transformation temperature, ease in 
purification, transparency, wider choice of materials, mini¬ 
mized convection effects and simplicity of experimenta¬ 
tion, organic systems (Glicksman et al 1983; Ecker et 
al 1989; Grugel and Poorman 1989; Podolinsky et al 
1989; Sangster 1994; Trivedi and Kurz 1994; Favite 
1996; Yasuda et al 1996) are better than metallic systems 
for scientific investigations. Succinonitrile is a low en¬ 
thalpy of fusion material and solidifies like a metal. In 
contrast, diphenyl is a high enthalpy of fusion material 
and solidifies like a nonmetal. As such a system involving 
succinonitrile (SCN) and diphenyl (DP) is an organic 
analogue of metal-nonmetal system giving eutectic and 
monotectic alloys. Preliminary investigations on this sys¬ 
tem have already been published earlier (Rai and Rai 
1996). However, in the absence of thermal data many 
aspects of this system could not be properly analysed. 
In view of this, phase diagram and thermochemistry of 
this system are reported here. 

2. Experimental 

2.1 Materials and purification 

Diphenyl obtained from Sojuz Chem., Moscow was 
purified by fractional crystallization with ethanol. Succino- 
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nitrile (Aldrich, Germany) was purified by repeated dis¬ 
tillation under vacuum. The purity of each compound 
was checked by determining melting point. 

2.2 Phase diagram 

Phase diagram of SCN-DP system was determined by 
the thaw-melt method (Rai and Mandal 1990). In this 
method mixtures of two components covering the entire 
range of composition are prepared in long-necked test 
tubes. The thawing and melting temperatures of each 
mixture were determined using a Toshniwal melting point 
apparatus attached with a thermometer which could read 
correctly up to ±0-5°C. A graph between melting tem¬ 
perature and composition gives the phase diagram. 

2.3 Enthalpy of fusion 

The values of enthalpy of fusion of each of the pure 
components, the eutectic and the monotectic, were 
determined by the DSC method using Mettler DSC~4000 
system. Indium metal was used to calibrate the system. 
Amount of sample and rate of heating were about 5 mg 
and 10°C/min, respectively, for each estimation. The 
maximum percentage error in the experimental value 
is 2. 

3. Results and discussion 

3.1 Phase diagram 

The phase diagram of SCN-DP system is given in figure 
1. The melting temperature of diphenyl decreases with 
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increase in concentration of succinonitrile. When con¬ 
centration of succinonitrile is 0-074 mole fraction, liquid 
immiscibility appears (Rai and Rai 1996). At this com¬ 
position which corresponds to 66-5°C melting temperature, 
two liquid phases, L, and L, and a solid phase are in 
equilibrium. This is the monotectic temperature and the 
composition corresponding to this temperature is the 
monotectic composition. When a liquid of monotectic 
composition is cooled below the monotectic temperature 
(T^), the monotectic reaction occurs and a liquid, L, 
which is rich in DP decomposes into a solid phase S 
(rich in DP) and another liquid L, (rich in SCN) according 
to the reaction, 

L, (rich in DP) ^ DP(5)-f (rich in SCN). (1) 

When composition of succinonitrile is increased beyond 
0-074 mole fraction, in a wide range area indicated by 
Lj -t- Lj, the two immiscible liquids are observed. At 
0-968 mole fraction of succinonitrile there is a eutectic 
point and at this temperature two solid phases and a 



Figure 1. Phase diagram of diphenyl-succinonitrilc system 
(o Melting temperature; • thaw temperature). 


liquid phase are in equilibrium. When a liquid of eutectic 
composition is allowed to cool below the eutectic tem¬ 
perature (Tg), the liquid decomposes to 

L, (rich in SCN) 6', (rich in DP) 

+ 5-2 (rich in SCN). (2) 

3.2 Thermochemistry 

It is well known that phase transformation involves two 
steps: (i) nucleation, and (ii) growth. While the nucleation 
mechanism depends on the solid-liquid interfacial energy, 
the growth mechanism depends on the manner in which 
particles from the liquid phase are added at the solid- 
liquid interface. The interface structure further depends 
on the entropy of fusion of the material under investigation 
and also upon the thermal environment in which the 
crystal is growing. Both the parameters, viz. interfacial 
energy and entropy of fusion can be calculated from the 
enthalpy of fusion data. Apart from this, radius of the 
critical nucleus, enthalpy of mixing, Jacksons roughness 
parameter and excess thermodynamic functions can also 
be calculated to throw light on the mechanism of 
solidification and formation of different microstructures 
of the eutectic and the monotectic. 

3.2a Enthalpy of mixing: The values of enthalpy of 
fusion of the pure components, the eutectic and the 
monotectic, determined by the DSC method are given 
in table 1. For the purpose of comparison the enthalpy 
of fusion of the eutectic, calculated by the mixture law 
(Rai and George 1996), is also given in the same table. 
If the eutectic is a simple mechanical mixture involving 
no interaction, the value of enthalpy of fusion determined 
experimentally and calculated by the mixture law should 
be the same. The enthalpy of mixing (A,^/f) (Rai and 
Shekhar 1994), which is the difference between the 
experimental value of the enthalpy of fusion and the 
calculated value of the enthalpy of fusion, is 0-4 kJ moT' 
in the present system. Thcrmochemical studies (Singh 
et al 1985; Rai et al 1987) suggest, that the structure 
of eutectic melt depends on the sign and the magnitude 
of the enthalpy of mixing. Three types of structures are 
suggested, quasi-eutectic for A^//>0, clustering of mole- 


Table 1. Heat of fusion, entropy of fusion and roughness parameter. 


Material 


Heat of 
fusion 
(kJ - mor') 


Entropy of Roughness 

fusion parameter 

(J • mol"') («) 
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Table 2. Interfacial energy of diphenyl, succinonitrile and their 
eutectic and monotectic. 


Parameter 

Value (ergem"^) 

SL, (DP) 

30-3 

SL, (SCN) 

9-3 

E (DP-SCN) 

10-0 

(DP-SCN) 

6-0 


cules for A^//<0 and molecular solution for A^//=0. 
The value of A^// being +04kJmor‘, suggests quasi 
eutectic structure (Rai et al 1987) in the eutectic melt. 

3.2b Entropy of fusion and excess thermodynamic func¬ 
tions: Theoretical studies on the entropy of fusion and 
calculation of some excess thermodynamic functions of 
eutectic predict the structure, stability and ordering in 
the eutectic melt. The values of entropy of fusion 
(Af^) of the pure components and eutectic were calculated 
using the equation, 

A,.// 

Af5=-^, (3) 

where A^// is the enthalpy of fusion, T the fusion 
temperature. The value of entropy of fusion being positive 
suggests that there is an increase in randomness during 
melting. Entropy factor is more effective in the eutectic 
in comparison to that of succinonitrile. On the other 
hand, it is less effective in comparison to that of diphenyl. 

The deviation from ideal behaviour and quantitative 
measure of molecular interactions can be expressed in 
terms of excess thermodynamic functions. They are given 
by the difference between the thermodynamic function 
of mixing for a real system and the corresponding value 
for an ideal system at the same temperature and pressure. 


Table 3. Radius of critical nucleus at different degrees of 
undercoolings. ‘ 


Undercooling 

(AT°C) 


Critical radius (xlO’^cm) 

DP 

SCN 

Eutectic 

2-5 


6-65 


3-0 


5-55 

4-62 

3-5 


4-75 

3-96 

4-0 


4-16 


4-5 


3-70 

3-08 

5-0 

2-25 


2-77 

5-5 

2-05 



60 

1-88 



70 

1-61 



90 

1-25 




= 


2oT 

m 

A^H-AT ' 


(5) 


where A^H and AT are melting temperature, heat 
of fusion and degree of undercooling, respectively. The 
interfacial energy is given by the relation 


CA.H 

a= -^- , ( 6 ) 

(yV)>^3(^jV3 

where A/ is the Avogadro number, the molar volume 
and parameter C lies between 0-30 and 0 33. The values 
of interfacial energy calculated using (6) are given in 
table 2. The critical size of the nucleus for pure com¬ 
ponents and the eutectics is given in table 3. It is evident 
from the reported data that for a sample, the size of 
the critical nucleus decreases with an increase in the 
undercooling. This is due to an increase in the amplitude 
of vibration resulting from an increase in the undercooling. 
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F(real)-A(ideal), (4) 

where Y is any of the thermodynamic functions. These 
functions can be calculated using the method reported 
earlier (Rai et al 1983). The value of being positive 
(210-7 J mol”') suggests (Wisniak and Tamir 1978) that 
associative interaction between like molecules is stronger 
than that between unlike molecules. 

3.2c Interfacial energy and size of critical nucleus: 
When a liquid is cooled below its equilibrium melting 
temperature, the liquid phase does not solidify sponta¬ 
neously because, under equilibrium condition, it contains 
clusters of molecules. So long as the clusters are all 
below the critical size (Christian 1965), they cannot 
grow to form crystals and so no solid is formed. The 
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Abstract. Depolarization current characteristics of solution grown pure ethyl cellulose (EC) films of about 
20 thickness have been studied as a function of electrode materials at constant poling field (5 x lO'* V/cm) 
and poling temperature 40°C. Thermally stimulated current (TSC) thermograms of EC consists of two well 
resolved peaks (located at 60°C and 140°C) for Al-Al system, which are attributed to the deorientation of 
strongly attached ethoxy groups of glycosidal units and diffusion of space charges either at electrodes or 
due to their thermal release at higher temperatures from the defect levels. For dissimilar electrode combinations 
(Al-Ag/Cu/Au/Sn/Pb), an indication of peak of lower magnitude at around (50-70°C) alongwith a higher 
temperature peak (140-155° C), have been observed. TSC parameters are found to change with the choice 
of electrode material. The dependence of dark current at 40°C in metal-ethyl cellulose-metal systems on 
applied voltage in the range (2 0-5 0) x 10“ V/cm has also been studied. The results of current-voltage 
measurement on EC have been interpreted to show that the Schottky-Richardson mechanism is the controlling 
transport mechanism. Zero field current density extrapolated from plots are found to vary with metal 

work function. 

Keywords. Thermally stimulated current; electrical conduction; ethyl cellulose; space charge relaxation; 
dipolar orientation; Schottky-Richardson mechanism. 


1. Introduction 

Thermally stimulated discharge (TSD) of polymers has 
become a widely used experimental technique for the 
investigation of various material parameters such as 
charge storage properties, determination of the mean 
depth of the internal charge, activation energies of traps 
and the trap structure of the material. The importance 
of carrier injection (Gutman and Lyons 1967; Koton 
1968; Lamport and Mark 1970) for the explanation of 
electrical properties of organic solids has been long 
recognized. It is only in recent years that the possibility 
of technological application of electrification of polymers 
due to their contacts with metal has been discovered 
(Shrivastava et al 1981). The availability of new mea¬ 
suring techniques such as photocurrent (Wintle 1974) 
and thermally stimulated discharge (TSD) current 
(Turnhout 1975) has freshened interest in achieving a 
better understanding and application of the phenomena. 
Many studies and reviews (Pulfrey 1972; Kryozewski 
1975; Khare and Srivastava 1992a, 1994; Khare et al 
1994a) have been devoted to the problem of contacts 
and surface charge influence but there is controversy 
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about the causes of the polarity and amount of charge 
transferred under given conditions. 

Ethyl cellulose (EC) is a thermally stable polymer that 
exhibits excellent chemical resistance and good mechani¬ 
cal properties and probably there is no report on contact 
electrification of EC using different metal contacts. Hence, 
it is worthwhile to investigate electret forming chara¬ 
cteristics of EC in contact with metals providing a range 
of work functions. In view of this an analysis of the 
thermally stimulated current (TSC) and steady state 
current-voltage (/-V) characteristics of EC polarized 
under identical conditions with a special reference to 
the nature of the electrode metals has been made to 
elucidate the electrode effect on injection in the polymer. 

2. Experimental 

The isothermal immersion technique was utilized for 
preparing thin films of ethyl cellulose (EC). The solution 
was prepared in a glass beaker by first dissolving 21 g 
EC in 30 ml of chemically pure chloroform at room 
temperature and continuously stirred for 30 min by means 
of a teflon-coated magnetic stirrer. Thereafter, it was 
stirred and heated to 40°C to yield a homogeneous 
solution. The glass beaker containing the solution was 


Ultrasonically cleaned vacuum metallized microscopic 
glass slides were immersed vertically into the solution 
for a period of about 30 min. After deposition of the 
film the glass slide was taken out and dried in an oven 
at 40°C for 24 h. This was followed by room temperature 
outgassing at 10"®torr for a further period of 24 h. 
Aluminium, copper, silver, lead, tin and gold metals 
were vacuum deposited on the upper surface of the film 
to obtain a sandwich configuration. The following elec¬ 
trode combinations were used—^Al-Al, Al-Ag, Al-Cu, Al- 
Pb, Al-Sn and Al-Au. The preparation of electrodes and 
measurement of thickness were the same as reported earlier 
(Khare and Srivastava 1992; Khare et al 1994b). The 
voltage was applied from a high voltage unit EC-4800 D 
and current was recorded with an electrometer (Keith ley 
600B), which was carefully shielded and grounded to avoid 
ground loops or extraneous electrical noise, at a linear 
heating rate of 4 K/min. The thickness of the sample was 
of the order of 20//m which was estimated by measuring 
the capacitance of the fabricated sandwiches taking the 
value of dielectric constant e of ethyl cellulose as 3. 

In the present case EC thermoelectrets were formed 
by a thermal method under identical polarizing conditions 
of temperature (7^, higher field and time (z^), viz. 
40°C, 5-OxlO‘^V/cm and 120 min but with different 
metal electrodes, namely copper, silver, lead, tin and 
gold, having different work functions. Thermoelectrets 
were shorted for 5 min to eliminate stray charges. The 
method of polarization of sample and measurement of 
TSC and steady current were same as reported earlier 
(Khare and Srivastava 1992; Khare et al 1993). 

Below 40°C it was difficult to measure the small 
currents generated by voltage application in the measuring 
set up. The film was kept at 40°C for about 2 h and 
then a voltage was applied and the transient behaviour 
of the current was noted. Voltage application gave a 
current which was found to decrease to a much lower 
value in a matter of 2-5 min and then very slowly attain 
an almost steady value in a few hours. The same 
procedure was followed to observe the steady value of 
current at each step of voltage increment. 

The depolarization kinetic data, activation energy {E), 
charge released {Q) and relaxation time (r„) for the 
observed peaks have been calculated using initial rise 
method of Garlick and Gibson (1948), Simson’s rule 
and Bucci et al (1966) theory, respectively (table 1). 

3. Results and discussion 

Thermally simulated current (TSC) thermograms of ethyl 
cellulose (EC) samples formed by incorporating alu¬ 
minium, copper, silver, lead, tin and gold as electrode 
metals with polarizing field, = 50 x 10“^ V/cm at 
T =40°C are shown in figure 1. In all the cases of 
electrode variation, an indication of peak at lower 


lying in between 14U-10D u nave been observed, ror 
Al-Al system, TSC thermograms are characterized with 
two peaks located at 60° and 140°C, respectively. Cur- 
rent-voltage (/-VO characteristics of EC at 40°C with 
polarizing fields (2-0-5-0) x 10"^ V/cm using different 
metal electrodes in I-E''- form are shown in figure 2, 
It is clear from figures 1 and 2 that the current through 
the sample is greater when the upper electrode aluminium 
is replaced by silver, copper, tin, lead or gold. 

The TSC thermograms are characterized by various 
peaks, due to dissipation of charges arising because of 
polarizations during charging process. The polarization 
of the material may arise due to various mechanisms, 
the important of which are dipole orientation, charge 
displacement, molecular and domain structures, macro/ 
microscopic displacement of ions with subsequent trap¬ 
ping, surface and space charge polarization etc. The 
trapping centres originate in the transition of the crys¬ 
talline state into amorphous regions, in the rearrangement 
of molecules in the polymer chain and from the presence 
of impurities and defects in the materials. The various 
defects, in addition to giving rise to their own chara¬ 
cteristic trapping levels, can also combine with other 
trapping centres present around them, and thus they can 
modify qualitatively as well as quantitatively the TSC 
spectra observed under polarizing conditions. EC is known 
to be a partial polar polymer. Its main chain is C-O-C. 
The side groups also have a structure of the form C-O-C 
but they differ in electronegative nature from the main 
chain. The electronegativity of the main chain bond 
arises from the difference in the electronegativity of 
carbon and oxygen atoms; however, the electronegativity 
of the side groups arises from the difference in electro¬ 
negativity of H 2 C- 0 -C 2 Hg and C-O-CjHj bonds in the 
side chains. The slight polar nature of the polymer is 
due to the difference in electronegativity of the main 
chain and the side groups. The dipolar relaxation is 
maximum around the glass transition temperature (T^). 
Its Tg lies in the temperature range 60-75°C for its 
atactic form and around 80°C for high molecular weight 
ethyl cellulose. In view of above characteristic of the 
polymer, it is expected to get first peak around 60-80°C. 
This peak seems to be due to the orientation of strongly 
attached ethoxy groups of glycosidal units. The relaxation 
process has been found to be associated with an activation 
energy of 0T8eV and related to the segmental rotations 
or translations of the main chain (Khare et al 1996). It 
remains unaffected with heating rate and storage time 
(Khare and Srivastava 1993). The peak current is almost 
linearly dependent on poling field. It is suggestive of 
the dipolar polarization (Khare et al 1992). The second 
peak appears on the higher temperature side of the TSC 
thermograms. This peak may be considered due to space 
charge polarization arising out of localization of bulk 
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generated and/or injected charge carriers in various availa¬ 
ble traps distributed energetically and spatially over the 
bulk of the sample. The joint motion/deorientation of 
the backbone chains of EC may be considered to con¬ 
tribute to some extent to this peak. The charge released 
from the trap levels at temperatures greater than are 
very well documented previously (Dubey et al 1990; 
Khare et al 1992; Khare and Srivastava 1993). An 
indication of peak of smaller magnitude (around 50-70°C) 
alongwith a broad TSC peak around 140-155°C have 
been found for dissimilar electrode combinations. These 
are space charge peaks modified by the type of metal 
(Shrivastava et al 1979), probably because the various 
metal-polymer interfaces possess different charge 
exchange rates, which change the space charge storage 
and current released by diffusion. The first stage of 
charging is carrier injection and the second stage is the 
entrapment of these charges in the border layer (Hino 
and Yamashita 1974). Consequently, the interpretation 
of depolarization appears to be reasonably possible in 
terms of charge detrapping of a space-charge build up 
due to carriers injected from the electrodes (Shatzkes 


1978). At high fields, the charges are injected from the 
electrodes into the polymer film and then trapped. It is 
only at temperatures above the glass transition of EC 
that the molecular chains are sufficiently agitated to 
release the charge stored in them (Khare et al 1994c). 
The amount and sign of charge injected depends on the 
relative work function of the metal-polymer interface, 
i.e. the difference in energy between the Fermi level in 
metal and the bottom of the conduction band in the 
insulator. Therefore, the TSC parameters are changed 
when the electrode forming material is changed. It also 
seems possible that the dipolar peak either gets masked 
by the space charge peak or its contribution to the total 
polarization is very small (Khare et al 1994a). 

The absence of peak at 60°C in Al-Sn and Al-Pb 
combinations may be due to oxidation of Pb and Sn 
metal electrodes. The ethoxy group being weakly elec¬ 
tronegative and the oxide layer electropositive, will neu¬ 
tralize each other’s effects to some extent. It is probably 
for this reason that in these two particular electrode 
combinations, the peak at 60°C is suppressed and only 
small humps are observed. 



Figure 1. Effect of electrode materials on thermally stimulated current for ethyl cellulose 
samples with polarizing field £p = 50x l(TV/cm at polarizing temperature Tp = 40°C. 
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In wide gap insulators and polymers with n-electrons, 
intrinsic carrier generation is not possible. Charge transfer 
from the metal depends on the electron levels in which 
the carriers shift freely under the influence of the applied 
field. The energy on this level may be given by (Baessler 
1972) 

= (4 +P'^ = 4-4 . (’) 

where is the ionization energy, the electron affinity, 
p* and p' are the polarization energies and ^ the ionization 
energy of a molecule of the dielectric. The observed 
relation between the current and voltage points to a 
conduction mechanism in which the charge carriers are 
released by thermal activation over a coulombic potential 
barrier that is decreased by the electric field. The potential 
barrier can be one of the two types; (i) barrier between 
electrodes and dielectric taking the classic image force 
into consideration (Schottky emission), and (ii) barriers 
due to trapping centres in the dielectric (Poole-Frenkel 
effect). Taylor and Lewis (1971) assumed a more genera¬ 
lized form of the potential barrier rather than a Coulombic 
barrier usually employed in treatment of the RS effect. 


They concluded that the potential barrier chosen is referred 
to the cathode-dielectric interface and probably is 
determined by a space-charge layer in the dielectric. 
Such space-charge layers are charged up by the absorption 
current. In the present case, plots may be interpreted 
in terms of Richardson-Schottky (RS) field assisted 
thermionic injection of carriers from metal electrodes. 

The primary source of charge carriers is the electrode 
contacts. Carriers are generated in the bulk on account 
of thermal ionization of impurities and absorbed water 
molecules. Nath and Kumar (1980) have shown injection 
of holes in cellulose acetate film. Hole injection results 
in depletion of layer of negative charges creating a field 
at the interface between the electrode and film surface. 
When two metals are brought in close contact, a potential 
difference of work function is created between 

the two metals. The semi-insulating polymer differs from 
a metal, however, in that an electric field may exist 
within the interior of the semi insulator. For this reason 
the contact potential drop between the metal and the 
semi-insulator may take place within the material rather 
than at the contact interfaces. Along with the field, there 
may exist a depletion in the accumulation of charges in 
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the surface layers (Khare and Srivastava 1992b). The 
work function of the metals used for contacts is lower 
compared to that of the work function of the polymeric 
material which creates an accumulation layer at the 
interfaces and leaves charged donors in the bulk. This 
in turn creates fields at the interfaces of the metal-polymer 
system (Khare et al 1994d). Suppose an insulator has 
a large energy gap with its Fermi level located above 
the mid gap and there is no surface state, then after 
contact with a metal surface, the electron transfer at 
contact equalizes both Fermi levels. The barrier height 
depends on the relative values of metal work function 
and the electron affinity of the insulator (Takai et al 
1975). The surface state (Wintle 1977) is one of the 


factors which affects injection phenomena. Donor-like 
surface states located between the conduction band and 
Fermi level, after contact tend to bend the band edge 
towards the Fermi level, making electron injection from 
the metal much easier. On the other hand, acceptor like 
surface states located between the valence band and 
Fermi level act as a stepping-stone for electrons injected 
from the valence band of the insulator to the metal, 
thus, enhancing hole injection (Shrivastava et al 1979). 
Electrons must overcome the potential barrier (x) given 
by 

X" = (p-A^, 

while holes are given by 



Table 1. Depolarization kinetic data of electrode variation in ethyl cellulose films 
(E = 5-Ox 10'^volts/cm, 7 =40°C). 

p p ' 


Electrode 

combination 

Peak 

temperature 

(°C) 

Peak 

current 

(A) 

Activation 

energy 

(eV) 

Relaxation 

time 

(sec) 

Charge 
released Q 
(Coul) 

Al-Al 


3-2 X 10-“ 

0-18 

4-789 X 10-“ 

6-222 X 10-“ 



1-5x10-"’ 

0-32 

8-698 X 10-"’ 

1-576 X 10-"’ 

Al-Ag 


7-2 X 10-” 

0-48 

9-978 X 10-"’ 

3-321 X 10"” 

Al-Cu 

145 

2-0x10-” 

0-52 

3-957x 10-” 

2-234x10'” 

A 1 A .. 

1 <rv 

'2.ri vy 1 n-8 

n.KA. 

<.000 vy 1 rv-i> 

Q.nAQ yy 1 0-8 
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= ( 2 ) 

where ^ is the metal work function, the electron 
affinity and ^ the ionization energy of a molecule of 
the dielectric. 

In the present case the magnitude of TSC and steady 
state current was found to be higher in dissimilar electrode 
than similar electrode combinations. This shows the effect 
of materials of electrodes on the current of the sample 
sandwiched between them. The observations are consistent 
with the proposed R-S mechanism in which the carriers 
are injected over the field-dependent polymer-electrode 
interfacial barrier. The emission current is described by 
Richardson’s equation 


= exp 


-1 

kT 


(3) 


V ; 

where is the current density, T the absolute temperature, 
k the Boltzmann’s constant and A a constant. 

The classical RS effect predicts a cur rent-voltage 
relationship of the form 


J~AT^ exp 


(- 

kT 

V y 


exp , 


(4) 


with 


^» kT 


4zc e € q d 


(5) 


where d is the film thickness, e the dielectric constant, 
the permittivity of the free space and e the electronic 
charge. For zero field, (4) reduces to (3) and the current 
density depends on the potential barrier at a constant 
temperature. The zero-field current density (J„) was extra¬ 
polated by extending the linear portion of the plot 
backwards to meet the current axis. The values of J„ 
vs metal work function have been plotted in figure 3. 
The negative slope of the line when EC makes contact 
with Ag which signifies that the electron injection 
decreases for Al-Ag electrode combination and so the 
zero field current decreases. The possible meaning that 
may be attached to the positively sloped line which is 
obtained when ethyl cellulose will be in contact with 
copper, tin, lead and gold is that, hole injection into the 
polymer film increases with the increase in metal function 
(Khare 1994) and hence zero field current also increases. 


4. Conclusion 

Thermally stimulated current and conduction current of 
ethyl cellulose films reveal a strong dependence on the 
electrode forming material which is the source of charge 
carriers responsible for conduction and observed polari¬ 
zation of the polymer. 
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Electrical properties of zirconium diselenide single crystals grown 
by iodine transport method 
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Abstract. Single crystals of zirconium diselenide (ZrSej) were grown by chemical vapour transport method 
using iodine as the transporting agent. The crystals were found to exhibit metallic behaviour in the 
temperature range 77-300 K and semiconducting nature in 300-443 K range. The measurements of thermo¬ 
electric power and conductivity enabled the determination of both carrier mobility and carrier concentration. 
The variation of carrier mobility and carrier concentration with temperature indicates the presence of deep 
trapping centres and their reduction with temperature in these crystals. 

Keywords. Zirconium diselenide; chemical vapour transport; electrical conductivity; thermoelectric power; 
trapping centres. 


1. Introduction 

Zirconium diselenide (ZrSe 2 ) crystallizes in a cadmium 
iodide type structure, which is characterized by strong 
bonding within layers and relatively weak Van der Waal’s 
bonding between adjacent sheets of layers. The structure 
is markedly anisotropic giving rise to an extended growth 
perpendicular to c-axis of the hexagonal crystals. These 
crystals have interesting electrical (Lee et al 1969a), and 
semiconducting (Lee et al 1969b) properties. They also 
exhibit some interesting switching and memory effects 
(Lee et al 1969a). A possibility of their use in solid 
state solar cell has also been predicted (Tributsch 1980). 

The data concerning electrical properties of ZrSe^ in 
the literature are very scarce and there have been no 
recent reports. It has been reported that ZrSe, grown by 
granular zirconium has a metallic nature whereas samples 
prepared from powdered zirconium show a semicon¬ 
ducting behaviour (Onuki et al 1982). They have also 
shown that in the temperature range 50-300 K the elec¬ 
trical resistivity of layered ZrSe, is proportional to T^. 
It has also been reported that samples with higher 
stoichiometry exhibit semiconducting behaviour while the 
less stoichiometric compounds are metallic (Deniard et 
al 1987). The normal semiconducting behaviour of ZrSe 2 
has also been reported (Zheng et al 1989) in electrical 
anisotropic studies in the temperature range 100-800 K. 
Looking at the anomaly in the resistivity behaviour of 
samples in the low and high temperature regions, we 
have carried out these measurements in both the ranges 


*Author for correspondence 


of temperatures. Moreover, the fact is ZrSe, is expected 
to have highly anisotropic electrical properties. We have 
also concentrated on this aspect in the resistivity mea¬ 
surements. In this paper, we, therefore, report the growth, 
characterization, measurements of d.c. resistivity, thermo¬ 
electric power, mobility and anisotropy in electrical 
resistivity in different temperature ranges. 

2. Experimental 

2.1 Growth of ZrSe^ single crystals 

Synthesis of zirconium diselenide was first reported by 
Van Arkel (1924). Later, Hahn and Ness (1957, 1959), 
McTaggari and Wadsley (1958), Green way and Nitsche 
(1965), Gleizes and Jeannin (1970), Remmington et al 
(1972), Lucovsky et al (1973) and Whitehouse et al 
(1973) have grown this compound by halogen transport 
methods. However, in the usual vapour growth experi¬ 
ments the quartz ampoule containing the charge is placed 
in the furnace in such a way that its bottom is directly 
in contact with the furnace wall while the upper portion 
is projected inside the furnace tube. This creates an 
uneven temperature gradient in the vertical direction 
leading to the introduction of defects into the grown 
crystals. In order to avoid this temperature gradient, 
arrangements were made to keep the quartz ampoule 
(along with the charge) coaxially within the furnace tube 
in the present investigation. The length of the quartz 
ampoule was about 250 mm with an inner diameter of 
20 mm. The hot and cold junctions of the furnace were 
maintained at 1123 K and 1073 K, respectively. Iodine 
was used as a transporting agent and large sized single 
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crystals of ZrSe, were obtained over a growth period 
of 500 h at the colder end of the ampoule. The crystals 
could easily be cleaved along the basal (0001) plane. 
The structural characterization and the compositional 
analysis were made by X-ray diffraction and energy 
dispersive analysis of X-rays (EDAX), respectively. 

2.2 Electrical resistivity measurements 


TEP (5'||) measurements were carried out in the temperature 
range 313-443 K. Due to experimental limitations, TEP 
(5'J measurements normal to the basal plane i.e. parallel 
to c-axis, could not be made. 

3. Experimental 

3.1 Structural characterization 


D.C. electrical resistivity measurements were performed 
on the freshly cleaved (0001) basal planes in a direction 
perpendicular to c-axis. Low temperature resistivity mea¬ 
surements along the basal planes Oo,,) in the temperature 
range 77-300 K were made in a specially designed 
cryostat LN-DP manufactured by Scientific Solutions, 
Mumbai. Electrical contacts were made using proper 
solder/silver paste. Four point contacts were used to 
obtain the resistivity. For resistivity (Pn) measurements 
in the temperature range 303-443 K a four-probe set up 
manufactured by Scientific Equipments, Roorkee was 
used. Resistivity measurements normal to the basal plane 
(pj were carried out on an experimental set up designed 
and prepared by the University Science and Instrumen¬ 
tation Centre (USIC), Sardar Patel University, Vallabh 
Vidyanagar. Starting from room temperature (303 K), the 
temperature of the sample was increased slowly in steps 
of 10 K until a temperature of 443 K was reached. At 
each step the corresponding value of resistivity was 
evaluated. 

2.3 Thermoelectric power measurements 

The thermoelectric power (TEP) measurements were ca¬ 
rried out using the thermopower (S) measurements set 
up TPSS-200 developed by Scientific Solutions, Mumbai. 
In all the measurements AT was kept as 4 K and the 


Figure 1 shows the zirconium diselenide (ZrSe^) single 
crystals obtained in the present investigation. Figure 2 
shows an a-axis oscillation photograph corresponding to 
the cleaved front face of ZrSe^ crystals taken with CuKa 
radiation employing a conventional X-ray tube with a 
camera radius of 3 cm and collimeter aperture of OT 
cm. The crystal was oscillated about the a-axis between 
25-40°. The a-axis oscillation photograph shows rows 
of constant h, k and varying 1. The presence of spots 
along the rows clearly indicates that the grown crystals 
have an ordered single crystal nature. The values of 
lattice parameters obtained from the X-ray diffraction 
are found to be a = 3-772 A and c = 6-122 A clearly 
matching with the reported values (Whitehouse et al 
1973) for the compound ZrSe,.yg+„. 3 . The composition 
of the crystals as obtained by EDAX correspond to 
ZrSe,y^. 

3.2 Variation of resistivity with temperature 

Results of low temperature (77-300 K) resistivity (/ 0 ||) 
measurements on a representative sample are shown in 
figure 3. It is seen that resistivity increases with increase 
in temperature and is found to be proportional to 7^ 
similar to TiSj and ZrSCj as reported by Onuki et al 
(1982). 

In the temperature range 303-443 K it is seen that 



Figure 1. Single crystals of ZrSe 2 . 






Figure 2. X-ray oscillation photograph of ZrSOj single crystal. 



Figure 3. Variation of resistivity with square of temperature. 



the resistivity decreases with an increase in temperature 
showing a semiconducting behaviour. The activation 
energies obtained from the graph (figure 4) are 0'0587 eV 
in the temperature range 303-393 K and 0-187 eV in the 
temperature range 393-443 K. 


3.3 Variation of electrical resistivity anisotropy with 
temperature 

Zheng et al (1989) measured anisotropy in the electrical 
resistivity (p^lpf), in the temperature range 100-300 K. 






In our investigation anisotropy measurements have been The electron density has been evaluated by using the 
carried out in the temperature range 303-443 K. It is above relation. The mobility, fi, of the charge carriers 

seen in figure 5 that the anisotropy ipJPy) ratio increases is determined from the relation 
with increase in temperature. This variation of anisotropy ^ 

with temperature shows an identical behaviour to that /< = —, (3) 

reported by Zheng et al (1989). 

where n is the carrier concentration, e the charge of 
electron and a the electrical conductivity. 

3.4 Variation of TEP with temperature 


The polarity of the thermally generated voltage at the 
hot end was positive indicating that the crystals are 
«-type in nature. The TEP was found to increase linearly 
with temperature. The TEP iS = S^) has been used to 
evaluate the carrier mobility using the relation given by 
Lee et al (1969b): 

TEP = 5 = -i{|.A), 
or 

5 = - {A + ln[2(2m* kT)^‘'^/nh^]}, (1) 

where A is the thermoelectric factor, n the electron 
density, h the Planck constant and m* the effective mass 
of electron. The above equation after the substitution of 
various constants and taking A = 1-5 for ionized impurity 
scattering simplifies to 

log n = (3/2) log r-0-0055'+15-719. (2) 



303 323 343 363 383 403 423 443 

Temperature (K) 

Figure 5. Variation of anisotropy with temperature. 


4. Discussion 

The dependence of electrical resistivity can be attri¬ 
buted to carrier scattering rather than carrier activation. 
However, due to experimental limitations the scattering 
and degeneracy could not be studied in detail in this 
range (77-300 K) of temperature. The high value of 
resistivity in the perpendicular direction suggests that 
conduction of electron in the direction normal to the 
layer is very little. Further, the anisotropic behaviour 
and its variation with temperature clearly shows that 
electronic conduction in ZrSe, is highly anisotropic. This 
large anisotropy in electronic behaviour may be attributed 
to the sufficiently larger effective mass anisotropy m* 
being 0-43 whereas m* being 7-00, where II and 1. 
denote, respectively, parallel and perpendicular to the 
layers of ZrSej single crystals (Isomaki et al 1979). The 
high values of resistivity and large anisotropy further 
suggest that electron states in ZrSe^ is two dimensional 
in contrast to a three dimensional behaviour seen in 
GaSe and SnSe. The low values of activation energy 
obtained from figure 4 can be attributed to the high 
carrier concentration density which gives rise to energy 
levels very close to the conduction band. These carriers 
will also be scattered by various centres or regions such 
as grain boundaries, dislocations, surface states etc in 
addition to the normally observed scattering modes, viz. 
ionized impurity, lattice or phonon, piezoelectric etc. The 
positive temperature dependence of mobility suggests 
that ionized impurity scattering might be the dominant 
scattering mechanism in these crystals in the observed 
temperature range. From the slope of temperature de¬ 
pendent mobility graph (figure 6), it is possible to obtain 
some qualitative idea about the scattering mode. Since 
all scattering processes are likely to follow the relation 
pi r*, the magnitude and sign of x will determine the 
scattering mode. In the present case the value of x from 
figure 6 comes to about 0-05 and is positive. The 
theoretical value of x for ionized impurity acting as a 
single mode of scattering is 3/2. The deviation in the 
observed value seems to be due to the superimposition 
of the ionized impurity scattering with some of the above 
mentioned scattering processes. The observed range of 
temperature i.e. 303-443 K, corresponds to the extrin¬ 
sic/impurity range only as seen from the low values of 
conductivity activation energy. Also, the carrier concen- 
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Figure 7. Variation of carrier concentration with reciprocal of 
tration can be expressed as temperature. 


n « exp 


(ii 

kT 


(4) 


V J 

where En is the carrier activation energy of the n-type 
ZrSej single crystal. The above equation predicts a linear 
relation between \ogn and 1/T as observed in figure 7. 
The carrier activation energy estimated from the slope, 
was found to be 0-29 eV. Since this value is much less 
than half the reported value of electrical band gap of 
1-22 eV (Murray et al 1972), the presence of deep donor 
like levels is possible in this material. Hence, heating 
in the lower temperature range reduces the number of 
deep trapping centres and thereby increasing the number 
of free carriers and carrier mobility. The mobility may 
reduce only at still higher temperatures where the sca¬ 
ttering due to lattice thermal vibration dominates. 


5. Conclusions 

Single crystals of zirconium diselenide could be grown 
by chemical vapour transport using iodine as the trans¬ 
porting agent. These crystals exhibited metallic behaviour 
in the temperature range 77-300 K and semiconducting 
nature in the temperature range 300-443 K. The positive 
temperature dependence of carrier mobility suggested 
that ionized impurity scattering might be the dominant 
scattering mechanism in these crystals. The variation of 
carrier mobility and carrier concentration with temperature 


in the semiconducting range indicates the presence of 
deep trapping centres. 
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Dielectric properties and electrical conduction in yttrium iron 
garnet (YIG) 
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Abstract. The dielectric properties (dielectric constant and loss) of a single crystal of yttrium iron garnet 
measured in the temperature range 77-725 K and in the frequency range 100 Hz-1 MHz. 
AC conductivity was derived from dielectric constant and loss. DC conductivity was measured in the 
temperature range 30-725 K. Thermoelectric power (TEP) was measured from 77-800 K. On the basis of 
the results, conduction in this garnet is interpreted as due to small polarons. The nature of conduction at 
different temperature ranges is discussed in the light of existing reports on defect formation. 

Keywords. Yttrium iron garnet; dielectric constant; conductivity; polaron conduction. 


1. Introduction 

Yttrium iron garnet (YIG) exhibits ferrimagnetism with 
low hysteresis loss and high resistivity at room tempera¬ 
ture. Because of this, it is of technical .importance and 
finds use in radio-electronics and UHF technology 
(Krupichka 1976). Hence the study of magnetic and 
electrical properties is important and useful. 

The magnetic properties of YIG have been investigated 
by many workers (Epstein and Frackiewicz 1959; 
Wurlitzer and Richter 1970; Wurlitzer 1987; Bush 1988; 
Chukalkin et al 1989). There are some reports on studies 
on electrical conductivity (Elwell and Dixon 1968; 
Fontana and Epstein 1971; Metselaar and Larsen 1974; 
Lai et al 1982; Petrov et al 1986). There are only a 
few reports on dielectric properties. Hirakata et al (1991) 
reported the dielectric behaviour at liquid nitrogen tem¬ 
peratures but the dielectric constant was not clearly 
indicated. Hofmeister and Campbell (1992) determined 
low- and high-frequency dielectric constants from infrared 
reflectivity measurements and indicated that the data by 
direct measurement is not available for a pure YIG 
single crystal. 

From resistivity and thermoelectric power measure¬ 
ments on HF-doped YIG in the range 500-1000 K, 
Elwell and Dixon (1968) concluded that the conduction 
should be described by localized model rather than the 
band model. Studies on electrical conduction in silicon- 
doped YIG by Fontana and Epstein (1971) in the tempe¬ 
rature range 20-200°C is interpreted as arising from 
localized hopping of electron amongst four Fe^"^ cation 
sites around Si'*'" impurity sites. 
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From DC conductivity and thermoelectric power mea¬ 
surements, Metselaar and Larsen (1974) indicate the 
possibility of p-type conduction in YIG in the high 
temperature range 1200-1700 K. Results of Lai et al 
(1982) and Petrov et al (1986) in high temperature range 
>600K and >1270K, respectively support conduction 
via thermally-activated hopping of p-type charge carriers. 

As evident from literature survey, while dielectric 
properties have been measured only in the liquid nitrogen 
range of temperature, resistivity measurements have been 
done on a Hf^’^'-doped sample at higher temperature 
beyond 500 K. Electrical conductivity studies have been 
made up to 475 K on a polycrystalline sample. The 
authors have neither come across a detailed investigation 
of dielectric properties over wide range of temperatures 
and frequencies nor conductivity measurements on a 
single crystal sample. 

It is well-known that the measurement of dielectric 
loss is a powerful method to detect relaxation due to 
formation of dipoles. AC conductivity in a wide range 
of frequencies combined with DC conductivity and 
thermoelectric power measurements helps to understand 
the conduction mechanism and identify the charge carriers 
responsible for conduction. 

Based on the above-mentioned facts, the measurements 
of dielectric constant and loss, DC conductivity and 
thermoelectric power have been carried out on a single 
crystal of YIG. The results are discussed to understand 
the conduction mechanism. 

2. Experimental 

All the measurements were carried out on a flux-grown 
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single crystal of YIG with parallel (111) faces. The 
sample was well polished and was in the form of a 
semi-circular disk of dimensions 0-67 cm^ x 0 0615 cm. 

In the low-temperature region, capacitance and loss 
were recorded using a fully automated computer- 
controlled CGA-83 transformer ratio arm bridge. The 
bridge was capable of measurement in three-terminal 
configuration. The accuracy of capacitance and loss is 
01% and 1%, respectively. The measurements were 
carried out at 17 frequencies in the range 10-100 kHz 
at temperatures from 77 K to 340 K with an interval of 
5 K and with an accuracy of 0-1 K. 

The high-temperature dielectric measurements were 
carried out using a Hewlett-Packard 4192 A impedance 
analyzer in conjunction with a laboratory-designed cell, 
in two-terminal method. The sample was coated with 
silver paint to ensure good electrical contact. Low- 
temperature resistivity measurements were carried out 
using an automated temperature controller and a Keithley 
digital electrometer. A Keithley electrometer (model 
610 C) was used for high-temperature DC conductivity 
measurements with an applied voltage of 3 volts. 

For low-temperature resistivity and thermoelectric 
power measurements, laboratory-built sample holders 
were used with a closed cycle helium cryo-refrigerator 
as the cooling device. The temperature gradient was 
controlled by using a Lakeshore temperature controller 
with an accuracy in temperature of 01 K. High- 
temperature thermoelectric power measurements were 
carried out using differential method and a Keithley 
nano-voltmeter for recording the data. 

3. Results 

3.1 Dielectric properties and conductivity 

3.1a Low-temperature measurements: The low-tempe¬ 
rature data was recorded as capacitance C{a)) and 
G(co), ctj as a function of temperature, where G((o) is 



Temperature (K) 

_^lkHz _*_5kHz ^lOkHz 30kHz __100kHz 

Figure 1. Variation of capacitance with temperature at different 
frequencies. 


conductance. The variations are to be treated as relative 
variations in dielectric constant and loss. The temperature 
variation of capacitance and loss at a few selected 
frequencies are shown in figures 1 and 2, respectively. 
The change in capacitance with frequency from 1-100 kHz 
is only 1% in low-temperature range and 2% around 
340 K. The variation with temperature is slow with 
distinct peaks around 130 K, 160K and a sharp increase 
at 280 K. This increase may be attributed to a large 
increase in ionic contribution to dielectric constant. 
Hirakata et al (1991) observed some anomalies in the 
real part of the dielectric constant and loss around 130 
K. They believe that these changes reflect the changes 
observed in the first-order magnetoelectric effect experi¬ 
ments of Kita et al (1988). Kita et al (1988) attributed 
the changes to some sort of phase transition, leading to 
a non-centrosymmetric structure below 125 K. 

Such a transformation could not be detected by X-ray 
and Mossbauer studies (Hirakata et al 1991). According 
to Hirakata et al (1991), the lattice strain due to the 
structural transition is so small that it could not be 
detected by X-ray, electron diffraction, Mossbauer and 
dielectric measurements. However, they support the 
existence of some sort of phase transformation, on the 
basis of the fact that the second-order magnetoelectric 
effect also shows a sharp decrease around 125 K (Takano 
et al 1991). 

The variation of loss (G/co) with temperature is shown 
in figure 2. The peaks around 130 K, 160K and 285 K 
are the manifestation of the effects observed in magnetic 
dispersion curves due to the presence of Si'*'^ in YIG 
(Wurlitzer and Richter 1970). Other studies on YIG such 
as magnetoelectric effect (Hirakata et al 1991) mentioned 
earlier and studies on magnetic effect by Torres et al 
(1993) indicate peaks around 130 K and 180 K. The 
effect is believed to be due to the electron transfer 
between Fe^"^ and Fe^^ ions creating local order and 
induced anisotropy. The temperature at which large in¬ 
crease is observed i.e. 295 K corresponds to 0/2 where 



_^lKHz _,_5KHz ^lOKHz _^.50KHz _^100KHz 

Figure 2. Variation of loss (G/oj) with temperature at different 
frequencies. 
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0 (600 K) is Debye temperature (Landolt-Bornstein 1991). of frequency at room temperature. It is well established 
AC conductivity is obtained from the data on that for ionic crystals an increase in s and loss (tan^) 
dielectric constant (s) and loss, tan (5 using the relation at low frequencies is observed for samples due to voids, 

dislocations and other defects. These spurious effects 
^Ac~^ ^ 0 ^ (1) cease to influence the value beyond a certain frequency. 

The frequency-independent value is taken as static 
where is the vacuum permittivity and oj the angular dielectric constant. For the sample of YIG used for 
frequency. Figure 3 shows plots of AC conductivity measurements, the dielectric constant is almost frequency- 
against reciprocal temperature. For clarity, the regions independent beyond 50 kHz at room temperature. Values 
of peaks in dielectric loss are not shown. It can be seen of the dielectric constant obtained in the present work 
that conductivity is frequency-dependent and tempera- and those reported by earlier workers are given in table 
ture-independent up to 250 K. The low-temperature DC 1. The values obtained from direct measurement are 
conductivity (figure 4) is found to be independent of close to one another. 

temperature. (Relative variation with temperature is to The variation of dielectric constant with temperature 
be noted and not the absolute value). at different frequencies is shown in figure 5. It can be 

In the magnetic loss studies in YIG, Wurlitzer and seen that large variations with temperature are observed 
Richter (1970) and Wurlitzer (1987) observed magnetic at low frequencies. The variation of loss with temperature 
loss dispersion when a small amount of silicon is present at different frequencies is shown in figure 6. The curves 
as impurity which gets incorporated during growth or show well-defined peaks. The peaks show a shift to 
as dopant. This is attributed to the electron hopping higher frequencies for higher temperatures indicating 
between Fe^^ and Fe^"" ions. The crystal used for the that they are relaxation peaks. The activation energy for 
present measurements is grown by flux in which divalent the relaxation process is obtained from the plot of log 
impurities Pb, Ca, Fe^^ and F are present in addition to peak frequency (ft>p) against reciprocal temperature shown 
Si'*^ (Haussuhl 1976). in figure 7. The activation energy is calculated to be 

1-01 eV. 

3.1b High temperature measurements: Initially, the di- The temperature variation of AC conductivity is shown 
electric constant and loss were measured as functions in the form of a log a vs T~' plot (figure 8). It can be 



Figure 3. Variation of AC conductivity with reciprocal tern- Figure 4. Variation of DC conductivity with reciprocal tem¬ 
perature at different frequencies. perature. 


Table 1. Dielectric constant of YIG. 



Frequency 

Dielectric 


Sample 

(temp.) 

constant 

Reference 

Polycrystalline 

2 MHz (R.T.) 

12 

Anderson (1959) 
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seen that a frequency-dependent slow increase is observed 
up to 500 K. The conductivity is found to be high for 
higher frequencies as in the case of variation at low 
temperatures discussed earlier. At 500 K, all the curves 
merge at a point showing frequency independence. The 
frequency dependence is shown once again beyond 500 K 
which is rather an unusual behaviour. These measurements 
were repeated to check for consistency. The same type 
of variation was observed in repeated measurements. The 
increase with temperature is faster beyond 500 K with 
higher values at lower frequencies and beyond 670 K 
an abrupt increase of more than two orders is observed. 

The DC conductivity in the temperature range from 
room temperature to 700 K is shown in figure 9. Room 
temperature conductivity is of the order of 10“'^ mho 
cm~^ It reaches an order of 10“'‘mhocm“‘ around 670 K 
and still higher beyond that temperature. It is to be 
noted that distinct kinks are observed at 500 K and 
575 K. The temperature at which the first kink is observed 
corresponds to the point at which a frequency-independent 
conductivity is observed in AC conductivity just at one 
fixed temperature. Perhaps the nature of conduction 



liX) 150 2i» 250 300 350 -100 

Temperature (‘’CJ 


_».lMHi 100 kHz 50 kHz 
_a_10kHz ^IkHz 500 Hz 

Figure 5. Variation of dielectric constant with temperature ai 
different frequencies. 



n 10 l:Hz __ 1 l:Hz 500 Hz 

Figure 6. Variation of tan <5 with temperature at different 


changes beyond this point as discussed earlier. The 
second kink at higher temperature corresponds to Neel 
temperature of YIG (573 K). Lai et al (1982) observed 
similar variation for several iron garnets at the Neel 
temperature. For a related material, DylG, Yadav and 
Lai (1979) have reported a break in the logo vs T"’ 
curve at magnetic ordering temperature. Change in the 
activation energy for conduction near curie temperature 
has been observed for Ti-substituted YIG by Om Prakash 
and Bahadur (1982). 

The significance of such an observation can not be 
explained at this stage. It can only be said that magnetic 
ordering has some effect on conductivity. In the entire 
region, conduction is extrinsic and AC conductivity is 
higher than DC conductivity. 

3.2 Thermoelectric power 

The variations of thermoelectric power with temperature 
in the low- and high-temperature ranges are shown in 
figures 10a and b, respectively. The thermoelectric 
voltage, A£' developed across the sample over a 



Figure 7. Variation of logrUp with reciprocal temperature. 



1 MHz __ 100 kHz __ 50 kHz ^ 10 !:Hz 
Figure 8. Variation of AC conductivity with reciprocal tern- 
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temperature gradient, AT is shown as S-AE/AT. 
Following the convention adopted by Lai et al (1982), 
the positive sign of S indicates negative charge carriers 
and vice versa. In the low-temperature range the thermo¬ 
electric power shows a gradual increase while decrease 
is observed beyond room temperature up to 700 K. 
Similar variation with an increase to a certain temperature 
has been observed in other oxide compounds (Om Prakash 
et al 1994). The charge carriers are n-type as shown 
by the sign of the thermoelectric power. The value 
remains almost constant at higher temperatures beyond 
700 K. 

4. Discussion 

YIG shows n-type conduction due to the presence of 
substitutional cations or anion vacancies. Doping of YIG 
with Si'*^ induces creation of Fe^"^ ions. The simultaneous 
presence of Fe^^ ions alter the electric and magnetic 
properties of garnets. The anisotropic magnetic resonance 
and slow relaxation phenomenon are examples of the 
electron exchange between the two irons (Judy 1966; 
Wurlitzer 1987). These effects in magnetic properties 
have been extensively studied as mentioned earlier. The 
sample used for the present investigation has marginal 
excess of donor concentrations. The sample charac¬ 
terization is given as Si‘‘'^ (0-003) + F (002) (Haussuhl 
1976) against the sum of divalent impurities Fe^^ 
(0-01) -f- Ca^"^ (0-001) •+• Pb (0-011) (concentrations are ex¬ 
pressed in the units of ions per formula unit of YIG). 
The presence of anion vacancies also introduces additional 
amounts of donor concentrations. 

YIG is a cubic crystal (I^3j-Oj‘’) with a structural 
chemical formula Y3Fe3'^Fe2‘*'0];. The trivalent iron cations 



Figure 9. Variation of DC conductivity with reciprocal tem- 
oerature. 


are distributed over two different sets of lattice sites; 
‘d’ sublattice in which each cation is tetrahedrally sur¬ 
rounded by oxygen ions and ‘a’ sublattice with octahedral 
coordination. The yttrium ions occupy ‘c’ sites charac¬ 
terized by a dodecahedral oxygen coordination. The Si'*'^ 
ions are known to structurally enter the ‘d’ sites forcing 
the reduction of some of the octahedrally-coordinated 
Fe^^ ions to Fe^^. The Fe""^ ion may be considered as 
Fe^^ cation plus an electron. The dipole formed with 
this electron and the positive charge located at tetrahedral 
‘d’ sites lead complex polarizability. This accounts for 
the relaxation phenomenon observed in dielectric loss. 
As the temperature is increased, some of the electrons 
get ionized and participate in electrical conduction. The 
dipolar component gives an additional term in AC field. 
The AC conductivity is found to be two orders more 
than the DC conductivity. Thus conduction is due to 
the localized charge carriers. Metselaar and Larsen (1974) 
have proposed a model for the band structure for YIG. 
They are of the opinion that in dealing with the band 




Figure 10. a-b. Variation of thermoelectric power with 
temperature. 
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structure of YIG the problems are similar to the case 
of transition metal oxides like NiO. Hence the conduction 
takes place in the narrow band derived from Fe (3^/) 
orbitals which are octahedral in character. Lai et al 
(1982) proposed energy band diagram for rare earth 
garnets. Fe^'^-Sd bands are shown to lie above the top 
of the uppermost completely filled 0^~-2p band. 

Under normal situation these garnets should exhibit 
high metallic conduction. Experimentally it is found that 
pure garnets are insulators with conductivity less than 
at room temperature. It may be concluded that 
Fe^"'-3d bands are localized. Hence in the presence of 
Si'*'^ impurity the mode of conduction could only be due 
to localized electrons. The conduction mechanism can 
be understood by treating it as narrow band semiconductor 
as in the case of transition metal oxides and rare earth 
compounds. 

Investigations of electrical conductivity in transition 
metal oxides (Adler and Feinleib 1970), rare earth 
sesquioxides (Subba Rao et al 1970; Lai et al 1976; 
Lai and Dar 1976) and in rare earth perovskites (Gschneid- 



Figure 11. Variation of log (cr^ - a^^,) with logw' al 90 K. 


ner and Eyring 1979) indicated polaron conduction 
through hopping mechanism. Though for DylG conduc¬ 
tion, thermally-activated hopping mechanism has been 
proposed by Yadav and Lai (1979). Om Prakash and 
Bahadur (1982) have explained the electrical properties 
of Ti-substituted YIG on the basis of electron hopping 
model as suggested by Austin and Mott (1969). For 
YIG, small polaron hopping mechanism is proposed by 
Petrov et al (1986) at high temperatures. 

Theoretical work by several workers over the years 
has provided some understanding of conduction in oxides 
and transition metal compounds (Bosman and Van Daal 
1970). For these materials, the interaction between elec¬ 
trons and optical phonons is strong and the conduction 
is explained on the basis of polarons. 

The treatment of conduction by polarons is discussed 
by several workers (Appel 1968; Austin and Mott 1969; 
Mott and Davis 1971). Polarons belong to two categories, 
large and small polarons. In the large polaron model, 
the conductivity is by band mechanism at all temperatures 
and the AC conductivity decreases with frequency. The 
small polarons conduct in band-like manner up to a 
certain temperature, the conductivity showing an increase 

-8 ,---, 





.> 4 .■; 6,7 .s y 1 

Loga>* 

Figure 12. Variation of log(a^-o-p^) with logw- at 200 K. 


Table 2. Polaron parameters. 


Crystal parameters 

Polaron parameters 

Calculated values 

(Static dielectric constant)“= 12-5 

Tp = 1/2 [.t:/6N]''^ 

0-918 A 

(High frequency dielectric constant)*’= 4-84 

cOq^Tji K0/h 

7-85X lO'^Hz 

Q (Debye temperature)'’ = 600 K 


7-900 

N (Number of sites per unit volume)'* = 160/a^ 

W =eV[8jreo£prp] 

0-993 eV 

Inter ionic distance^ = 5-7 A 

m* = 20 £■ /r CO /e^ 

P P 

5-88X lO^m 

e 

a (Lattice constant)* = 12-376 A 

a = e^ {m*l2h^ o)^ 

13-85 


//Ip = 0-02 m* 

4-33X lO^m 

C 


“Present value, '^Hofmeister and Campbell (1992), ‘'Landolt-Bornstein (1991), ‘'Geller (1960), * Austin 
and Mott (1969) and ‘Gschneidner and Eyring (1979). 
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with frequency. At higher temperatures, the conduction 
is by thermally-activated hopping mechanism. The 
experimental results in the present case also shows the 
trend expected for small polaron conduction mechanism. 
It has been shown (Appel 1968; Adler 1968) that for 
ionic solids the concept of small polaron conduction is 
valid. 

As an evidence for small polaron conduction, we have 
examined the variation of conductivity with frequency 
on the basis of relaxation process, as applied to bound 
small polaron (Adler and Feinleib 1970; Bosman and 
Van Daal 1970). Adler and Feinleib (1970) have shown 
that for conduction by small polarons the conduction 
increases with frequency and the following relation holds 
good: 


where co is the angular frequency and t the staying time 
(-10"“'sec) for frequencies a)^r~<\. log (a^^ - Oq^,) vs 
logft>^ should be a straight line. Plots of log (cr^ - (Joj.) 
as a function of co^ at two temperatures 90 K and 200 K, 
are shown in figures 11 and 12, respectively. It can be 
seen that the plots are straight lines. This indicates that 
the conduction is due to small polarons. 

The important parameters associated with polarons are: 
polaron radius, r^, polaron binding energy, W^, effective 
mass, wip and dimensionless electron-phonon coupling 
constant a. The polaron parameters were evaluated from 
the data on dielectric constant and other crystal parameters 
using the relations given by Austin and Mott (1969). 
The input data and the polaron parameters are shown 
in table 2. 

For small polarons r^ should be less than the inter 
ionic distance and the coupling constant should be greater 
than 6 (Bosman and Van Daal 1970). The values of the 
parameters are appropriate for small polarons. 

5. Conclusions 

The static dielectric constant for YIG is found to be 
12-5 at room temperature at 1 MHz frequency. The 
disturbances observed in the dielectric constant around 
130 K, 160 K and in the dielectric loss at three tem¬ 
peratures 130 K, 160 K and 285 K correspond to the 
peaks observed in magnetic loss and induced anisotropy 
measurements on silicon-doped samples. No indication 
of phase transition is observed. From AC and DC 
conductivity measurements, creation of Fe^'^ due to the 
presence of Si'*'" as impurity and electron hopping between 
Fe^'^ and Fe^^ ions is believed to be the cause for 
conduction. Thermoelectric power data shows n-type 
charge carriers up to 700 K. 


temperatures and by thermally-activated hopping mecha¬ 
nism at high temperatures. The relaxation pattern observed 
in the dielectric loss and the frequency variation of 
conductivity at low temperatures support this conclusion. 
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Optical and mechanical characteristics of zinc sulphide-thorium 
fluoride mixed composition thin films for use in the near infrared 
region (1-10 pim) 
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Abstract. Mixed composition thin films of zinc sulphide-thorium fluoride have been deposited on glass and 
silicon substrates by thermal evaporation of mixtures of these materials in different proportions, from a 
single resistively heated source. The films are characterized for their optical properties (refractive index and 
extinction coefficient), mechanical properties (intrinsic stress), surface morphology and chemical composition. 
It is found that these films have tailorable refractive indices and low losses, and that fllms with certain 
compositions have low intrinsic stress and smooth surface morphology, making them suitable for incorporation 
in thin film multilayers for use in the near infrared region up to at least 10/im. 

Keywords. Mixed composition thin films; zinc sulphide-thorium fluoride; optical and mechanical chara¬ 


cteristics; near IR region. 

1. Introduction 

It is well known that there are relatively few single 
materials available which are suitable for incorporation 
.in thin film multilayer stacks useful in the near infrared 
region (l-lO/^m). The fabrication of multilayer optical 
coatings with specified spectral characteristics in this 
region would be greatly simplified if low loss optical 
materials with tailorable refractive indices (especially in 
the low index range) and acceptable mechanical properties 
(low intrinsic stress, adequate hardness, etc) are available. 
Mixed composition films offer the possibility of solving 
this problem. So far, most studies of mixed composition 
films have been confined to the visible region. Only in 
recent years has interest been shown in mixed composition 
films useful in the near infrared region (Sankur et al 
1988, 1991; Targrove and Murphy 1990; Gluck et al 
1992, 1993). The mixed composition films reported so 
far, along with their regions of transparency, are listed 
in table 1. 

Zinc sulphide and thorium fluoride are high and low 
index materials which are used for fabrication of multi¬ 
layers in the near IR region. We therefore investigated 
mixed composition thin films of zinc sulphide and thorium 
fluoride in the search for a suitable thin film material 
with desirable properties. 

2. Deposition of thin films of ZnS-ThF^ 

Both zinc sulphide and thorium fluoride have the same 


evaporation temperature (1100°C). Therefore, mixtures 
of the two materials could be evaporated from the same 
resistively heated source. The two materials were taken 
in different proportions by weight (measured by a Mettler 
microbalance), ground and mixed together thoroughly in 
a mortar and pestle and put into a tungsten boat, covered 
with a perforated shield to prevent spattering of the 
material. The thermal evaporation was carried out under 
a vacuum of about 3xl0"'^mbar, and the thickness of 
the deposited film was monitored by an optical thickness 


Table 1. Mixed composition films reported in the literature. 


Composition 
of film 

Region of 
transparency (am) 

Reference 

ZnSe-LaFg 

0-8-2-0 

Gluck era/(1992, 1993) 

ZnSe-SrFj 

0-6-80 

Sankur et al (1991), 
Gluck et al (1993) 

Si-YFg 

T5-7-0 

Sankur et al (1991), 
Gluck et al (1993) 

Si-ZnSe 

1-5-7-0 

Gluck era/(1992, 1993) 

Si-ZnS 

1 •5-7-0 

Gluck et al (1993) 

LaFj-BaFj 

0-25-2-0 

Targrove and Murphy 
(1990) 

Ge-MgFj 

2-0^-0 

Sankur et al (1988) 

Ge—CeFj 

2-0-5-0 

Sankur et al (1988) 

Ge-ZnS 

2-0-12-0 

Sankur et al (1988) 

Ge-CdTe 

2-0-100 

Sankur et al (1988) 


*.1. 
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monitor to keep track of the film thickness. Films of 
physical thicknesses in the range 0-3-1-O/rm were de¬ 
posited on substrates of ophthalmic glass (for refractive 
index, extinction coefficient and thickness measurements, 
as also for SEM and XRD studies), thin glass microslides 
(for stress measurements) and polished Si wafer pieces 
(for IR transmittance measurements and XPS studies for 
the chemical composition). 

3. Determination of refractive index, extinction 
coefficient and thickness of the fdms 

The refractive index, extinction coefficient and thickness 
of the films were determined from the trace of the 
transmittance vs wavelength of the films by the well 
known procedure described by Swanepoel (1983). The 
transmittance of the ophthalmic glass substrate coated 
with the film was measured as a function of wavelength 
from about 400-1500 nm. Depending on the thickness 
of the film, a number of maxima and minima were 
observed in the spectrophotometric trace. Figures la and 
b show typical traces for a pure ZnS film and a mixed 
composition ZnS-ThF^ film deposited on glass substrates. 


ZnS film on glass (ZnS-1) 



The envelopes of the maxima and minima in these curves 
were drawn, and the values of transmittance (from actual 
transmittance curve or from envelope curve) at the wave¬ 
lengths corresponding to each maximum or minimum 
were noted. A computer programme was developed which 
calculated the refractive index, n, extinction coefficient, 
k and film thickness, d at each of these wavelengths, 
on the basis of these transmittance values. These calcu¬ 
lated results are also shown in figures la and b, along 
with the transmittance curves. The variation in the cal¬ 
culated values of d at the different wavelengths, and the 
corresponding variations in the refractive indices calcu¬ 
lated using these values of d, gave estimates of the 
uncertainties in the values of n and d at each wavelength. 
The values of k are of the order of 10"^ which shows 
that these films are almost transparent in the near infrared 
region. 

The refractive indices, extinction coefficients and thick¬ 
nesses of a few selected films were also determined 
from ellipsometric measurements at 546-1 nm wavelength, 
and were found to agree with the values determined by 
the above method to within a few per cent. The thick¬ 
nesses of a few films were also determined directly by 



Figure 1. Plots of the transmittance vs wavelength for two films: ZnS-l (a) and ZnSThF-4 (b), deposited on ophthalmic glass 
substrates. The envelopes of the transmittance maxima and minima are shown. The results of calculations (§ 3)—the values of 










a stylus scan measurement in a Talystep instrument, and 
the values were again found to be in good agreement 
with the previously determined values. 

In figure 2, we have plotted the dispersion curves for 
the refractive indices of a pure ZnS film and a film 
(ZnSThF-4) fabricated from a mixture with 60% ZnS 
and 40% ThF^ by weight. The refractive index shows 
a regular decrease with increasing wavelength in both 
cases, as expected. 

In figure 3, we have plotted the refractive index of 
each film at 1 /<m (obtained by interpolation from the 
refractive index vs wavelength data) against the ratio / 
for the film, defined as 

/=W/(1+H'), , 

where w^^ is the weight of ThF^ (low index, n^^, constituent) 
in evaporation mixture and the weight of ZnS (high 
index, njj, constituent) in evaporation mixture. 

As we shall show later (§ 7), the fractional composition 
by weight of a film is found to agree with the fractional 
composition by weight of the evaporation mixture. 

There are various models (linear and non-linear) that 
predict the refractive index of a mixed composition film 
(Niklasson et al 1981), and some researchers, e.g. Chen 
et al (1996), have tried to fit their refractive index vs 
fractional composition data for mixed composition films 
to one or the other model. However, the scatter in the 
refractive index vs / data in figure 3 does not warrant 
the fitting to a particular model. Other researchers (Gluck 
et al 1992) have also expressed the same opinion. We 
can, however, say that in the present case the refractive 
index of the ZnS-ThF^ mixed composition film varies 
monotonically with composition (represented by the frac¬ 
tion, f) and follows an approximately linear relation. In 
fact, the dotted line in figure 3 is the result of the 
Linear Model (Niklasson et al 1981), according to which 
the effective refractive index of the mixed composition 
film is given by 



WavelengLK (nm) 

Figure 2. Rclfaclive index (calculated in § 3) vs wavelength 
for the twh” films ZnS-1 and ZnSThF-4. 


/i,ff = (/Ih • Wh + • vVl)/(w„ + wj 

(1 f') ■ ’ 

using the expression for / given earlier. 

4. Stress measurement 

The intrinsic stresses in the mixed composition films of 
ZnS-ThF^ are expected to be small, since ZnS films 
have been reported to have moderate compressive stress 
and ThF^ films have moderate tensile stress. The intrinsic 
stresses in the films were determined by a laser scanning 
set-up. Light from a HeNe laser was passed through a 
beam expander and aperture to obtain a collimated beam 
which would not suffer much divergence over a distance 
of several metres. The laser beam was split by a beam¬ 
splitter into two beams: one was sent off to the side, 
and the other was directed to the sample under inves¬ 
tigation. After reflection from the sample, this beam was 
directed by the same beamsplitter to a position sensitive 
quadrant detector, or to a screen a few metres away. 
The laser beam spot on the sample was translated over 
the surface of the sample by a micrometer over a distance 
of several mm, in a direction normal to the incident 
light, and the resultant deflection of the reflected light 
beam (which would indicate curvature of the sample 
surface) was measured by the position sensitive detector 
or on the screen. For deflections of less than 1 cm, the 
detector could measure the deflection to an accuracy of 
01 mm. But greater accuracy was usually obtained by 
measuring a deflection of a few cm on the screen. 

The intrinsic stress in the film was determined by the 
following relation; 



O 0-5 1-0 

fure ZnS Weight Fraction. Tore ThF^ 

Figure 3, Refractive index at 1 //m (calculated in § 3 and 
interpolated to 1 ^tim) vs the weight fraction, for various mixed 
composition films. The dashed line is the plot corresponding 
to the theoretical Linear Model. 




Film stress a= ■ - ■■ —— 
6tfO 


1 

R 



where E is the Young’s modulus, the Poisson’s ratio, 
t the thickness of substrate, L the film thickness and R 

S 1 

and Rq are the radii of curvature of substrate before and 
after deposition of film, and 

R = 2xL/d, 

where x is the displacement of , sample, L the distance 
of screen/detector from sample and d the displacement 
of beam spot, and similarly for R^. 

For the thin glass microslide used as the substrate for 
stress measurements (Corning 0211 alkali zinc borosilicate 
sheet glass), we used the values [Corning catalog]: 
£; = 745x 10'"N/m^ /< = 0-22, t =0-15 mm. 

The radius of curvature of the substrate was measured 
at the centre of the substrate before and after deposition 
of the film. The results of measurement are shown in 
figure 4. Compressive stresses are shown as positive 
and tensile stresses as negative. For ease of representation 
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Figure 4. Intrinsic stress (in kg/cm^) vs the weight fraction 
for the various mixed composition films studied in § 4. Com¬ 
pressive stresses are shown as positive and tensile stresses as 
negative. 


of positive and negative stresses in the same figure, and 
since stresses have usually been expressed in kg/cm^ in 
the literature, we have expressed the stresses in our 
films in kg/cm^, using the conversion 1 kg/cm^ = 9-8 x 10“* 
N/m^. We see that as the weight fraction of ThF^ in 
the mixed composition films increases, the intrinsic stress 
changes from compressive to tensile, and a film with a 
weight fraction of about 0-5, i.e. a film with about 
50 : 50 ratio of ZnS-ThF^, has very low stress. This is 
an illustration of the reduction of intrinsic stresses in 
mixed composition films reported by several researchers 
(Austin 1973; Sankur et al 1988, 1991; Targrove and 
Murphy 1990; Gluck et al 1992, 1993). 

5. Surface morphology 

The mixed composition films were examined under a 
scanning electron microscope to see the surface roughness, 
cracks or other surface features of the films. The films 
were examined both at normal incidence and at oblique 
incidence up to a magnification of 15,000 x. However, 
all the films exhibited featureless surface morphology 
up to this magnification. We therefore anticipate that 
these films will exhibit very low optical scattering losses 
and hence are suitable for use in optical multilayer stacks 
without introducing appreciable scattering losses. 

6. Transparency in near IR (1-10/^m) range 

The IR transmittance spectrum of thin films coated on 
silicon substrates was measured on a Perkin-Elmer spec¬ 
trophotometer from 2-5-10/rm. The spectra showed 
absorption peaks of the silicon substrate at about 
9-1 fim, and were otherwise quite flat over this entire 
range. We therefore conclude that these thin films are 
almost transparent over the wavelength range of interest, 
and are thus suitable for the fabrication of multilayer 
stacks. 

7. Structure and chemical composition 

X-ray diffraction studies of some of the mixed compo¬ 
sition films showed a peak corresponding to zinc sulphide, 
but there was no defined peak corresponding to thorium 
fluoride. This is clearly seen in the plots in figure 5, 
for films ZnSThF-1 (a), ZnSThF-4 (b) and ZnSThF-8 
(c). The films in (a) and (b) have large weight fractions 
of ZnS, and the diffraction peaks corresponding to ZnS 
are clearly seen. The film in (c) is predominantly ThF^, 
and no diffraction peak is seen here. This indicates that 
zinc sulphide exists in a crystalline phase in the mixed 
composition films, while thorium fluoride exists in an 
amorphous phase. • " 

The chemical composition of the mixed composition 
films was determined by the XPS technique, where 
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Figure 5. X-ray diffraction plots for three mixed composition 
films, showing the peak corresponding to ZnS and the absence 
of a peak corresponding to ThF^. 

surface scans of some of the films showed that the films 
were comprised of ZnS and ThF^ only. The weight 
proportions of ZnS and ThF^ in the films were found 
to be quite close to the weight proportions in the 
evaporation mixture. 

8. Conclusions 

The mixed composition films of zinc sulphide and thorium 
fluoride have tailorable refractive indices (from about 
2-3-1-5 at I pm wavelength) and fairly low extinction 


coefficients, of the order of 10'^ or lower. The intrinsic 
stress is very low for film compositions with a 50:50 
ratio by weight of ZnS and ThF^. The films exhibit 
almost featureless surface morphology and therefore 
should have low optical scattering losses. Further, while 
ZnS appears to exist in a crystalline phase in these 
films, ThF^ appears to exist in an amorphous phase. In 
conclusion, these mixed composition films should be 
useful in the fabrication of multilayer stacks to meet 
specified spectral characteristics in the near IR region. 
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Abstract. Magnetization measurements of the ferrite system Mgj^^Mnj^Fe 2 _ 2^04 were undertaken primarily 
to study variation of the saturation magnetization with manganese concentration and to gain information 
about the type of magnetic ordering. The observations indicate the existence of Yafet-Kittel (YK) type of 
magnetic ordering in this system which is well supported by results of Mossbauer effect measurements. A 
molecular field analysis of the YK spin ordering using the three-sublattice model is shown to explain the 
experimental data satisfactorily. 

Keywords. Ferrites; magnetic ordering; collinear; non-collinear. 


1. Introduction 

This possibility of changing continuously the concentra¬ 
tion of non-magnetic ions in the different sublattices 
makes the ferrites very interesting materials to study the 
appearance of various magnetic structures (Bhargava and 
Zeeman 1980; Brand et al 1985; Muralidharan et al 
1985; Dormann et al 1987; Chakravarthy et al 1988; 
Dormann and Nogues 1990; Jotania et al 1992). Among 
these, systems with spinel structure seems to be particu¬ 
larly attractive, as they allow a variety of magnetic 
disorder and frustration to be introduced. This is due to 
the fact that, in spinel intra-sublattice interactions are 
weaker than the inter-sublattice interactions, as a result 
there are unsatisfied bonds in the ferrimagnetic phase. 
Because of these unsatisfied bonds increasing magnetic 
dilution accentuates the competition between the various 
exchange interactions resulting in a variety of magnetic 
structures (Bhargava and Zeeman 1980; Brand et al 
1985; Muralidharan et al 1985; Dormann et al 1987; 
Chakravarthy et al 1988; Dormann and Nogues 1990; 
Jotania et al 1992). 

A large number of investigators have studied Ti- 
substituted spinels in order to understand their magnetic 
properties. Brand et al (1985) found in the case of 
that on increasing dilution x, the collinear 
ferrimagnetic phase breaks down before reaching the ferri¬ 
magnetic percolation threshold and for higher concentration 
of Ti, system approaches toward spin-glass. Dormann et 
al (1987) have also studied the magnetization measurements 
of Ti"*^ substituted lithium ferrite using .different techniques 
and suggested the possibility of canted spin arrangements 
to explain the magnetization measurements. 
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In order to understand the influence of magnetic ion 
Mn'*^ in place of non-magnetic ion Ti'*^, we have modified 
the magnetic phase of Mg,^^Ti_^Fe 2 _ 2^04 by replacing 
Ti'*"' by Mu'*"”, so that intra-sublattice interaction 
will become predominant over the inter-sublattice inter¬ 
action (/ab), which may lead to a frustration in the spin 
sub systems. In our recent communication (Jani et al 
1997), we have reported magnetic properties on the 
mixed spinel Mgj^^Mn/^ej.j^O^. The present paper 
reports the results of X-ray, magnetization and Mossbauer 
effect measurements on the Mg,^^Mn_^Fe 2 ^ 2>:^4 system. 
The variation of the saturation magnetization per formula 
unit in Bohr magneton number (n^) with Mn-concentration 
in this system is similar in nature to that in the Ni-Zn 
(Satyamurthy et al 1969) and Cu-Zn (Kulkarni and Patil 
1982) ferrites and the existence of Y-K angles on the 
B-site, observed in the latter systems is strongly suspected. 
A molecular field analysis of the Y-K spin ordering 
using a three-sublattice model is shown to explain the 
experimental data satisfactorily. 

2. Experimental 

Nine samples of Mn-substituted Mgi^^Mn^Fej^^^O^ sys¬ 
tem were prepared by the usual double sintering ceramic 
method for 0-l<jr<0-9 in steps of 0-1. The starting 
materials were analytical reagent grade oxides, FejOj, 
MgO and MnOj. These oxides were mixed in proper 
proportions and presintered at 990°C for 24 h. In the 
final sintering process the material was held at 1050°C 
for 12 h and slowly cooled to room temperature (2°C/min). 
The X-ray diffractograms were obtained using FeK„ 
radiation on a Philips X-ray diffractometer model PW- 
1820. The diffraction patterns show the sharp lines 
corresponding to a single phase spinel for all the samples. 
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The magnetization measurements of each sample at 
300 K were carried out using high-field hysteresis loop 
technique (Radhakrishnamurthy et al 1971). The Moss- 
bauer spectra were obtained at 12 K and 300 K in 
transmission geometry with source of lOmCi ^CoCPd) 
and a constant acceleration transducer interface to PC- 
based 1024 multi channel analyser. 

3. Results and discussion 

All the compositions of the system Mg,^,.Mn^Fe 2 _ 2 ^P^ 
exhibit single phase cubic spinel structure. No reflections 
other than those belonging to a spinel structure were 
observed in the patterns (figure 1). The values of lattice 
constant ‘a’ (A) determined from X-ray data with an 
accuracy of ±0-002 A for jc = 0-1 to 0-9 are shown in 
figure 2 as a function of x. It is clear from figure 2 
that the ‘n’ parameter initially increases up to x = 0-4 
and thereafter it levels off for further increase in 
a: >0-5. The observed nonlinear behaviour in ‘a' with x 
may be attributed to the simultaneous replacement of 
Fe^^ (0-64 A), Mg-" (0-60 A) and Mn'" (0-60 A). The 
cation distribution of MgFe^O, (a: = 0) has been reported 
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Figure 1. Typical X-ray diffraciograms of the system 
Mgi+^Mn^Fe 2 _T ^04 for a = 0-2, 0-6, 0-8 and 0-9. 


to be 

(Mg„,«Fe„,„)nMg,„,Fe,,«l“0„ (1) 

where A is the tetrahedral site and B the octahedral 
site. In order to determine the cation distribution, XRD 
intensities were calculated using the formula suggested 
by Buerger (1960): 

= ( 2 ) 

where notations have their usual meanings. The distri¬ 
bution of divalent, trivalent and tetravalent cations 
amongst octahedral and tetrahedral sites in the 
Mg| ^^Mn^Fe 2 _ 2,04 samples are determined from the ratios 
of intensities of X-ray diffraction lines: ^ 221 /^ 40 ’ A 2 (/A()o 
and / 4 ()(// 422 - t^ble 1, the results of X-ray intensity 
calculations for typical samples a =0-2, 0-6, 0-8 and 0-9, 
are listed along with experimental intensity ratios. It is 
evident from table 1 that tetravalent Mn-ions occupy 
B-site replacing B-site iron and the added Mg'" ions 
substitute for A-site iron without any change with x in 
the initial degree of inversion i.e. (1). 

In deriving the cation distribution of Mg, ^^Mn^Fe 2 _ 2 P 4 
our X-ray intensity calculations (table 1) and relative 
Fe^" intensities of Mbssbauer spectra (table 2), for A 
and B site ions were taken into account. Accepting the 
cation distribution of MgFCjO^ (1), an approximate cation 
distribution for the Mg, ^^Mnpe 2 _ 2 ,P 4 can be written as 



Figure 2. Variation of lattice parameter a (A) with Mn con¬ 
centration X. 
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Table 1. Comparison of X-ray intensities ratio for Mg, ^^.Mn^FeT_2^^0^ in estimating cation distribution. 


X 

Lattice 
constant 
a (A) 

A-site cation 

B-site cation 

42(/^440 

42(/4(X) 

^40(1 

/^422 

■ Obs. 

Cal. 

Obs. 

Cal. 

Obs. 

Cal. 

0-2 

8-3853 

^S().2xE^1).72 

MSoyiMn^, 2()Fe(, 

0-82 

0-82 

1-33 

1-35 

2-49 

2-14 

0-6 

8-4320 

^§()fi6^®0-34 

Mg(,.y4Mn„ g(,Fej, 

0-78 

0-53 

0-76 

0-66 

6-10 

4-46 

0-8 

8-3984 

^So-82^®().|K 

j,„Fe(, 22 

0-74 

0-33 

0-66 

0-43 

6-80 

6-30 

0-9 

8-3965 

Mg(,.yiFe(,()y 

Mgo.yyMn(,y()Fe(,,,, 

0-81 

0-29 

0-74 

0-35 

6-23 

7-87 


*fl(A) error ± 0-002 A. 


Table 2. Mossbauer parameters of Mg, ^^.Mn^pe, obtained at I2K. Average sublattice magnetization, /j.ix) in 


X 

Isomer shifts (mm/s)* 

Hyperfine fields (KOe) 

Mossbauer 
intensity ratio 

/r(x) = «^(x) 
in 


/5b 



0-2 

0-37 ±0-05 

0-37 ±0-05 

506 ±3 

523 ± 3 

1-22 

1-39 

0-6 

0-32 ±0-05 

0-32 ±0-05 

489 ±3 

510±3 

1-35 

2-50 

0-8 

0-32 ±0-05 

0-32 ±0-05 

469 ±3 

495 ±3 

1-22 

2-98 

0-9 

0-32 ±0-05 

0-32 ± 0-05 

448 ±3 

480 ±3 

1-22 

3-15 


‘*‘15 w.r.t. «'Fe. 

(MS().08 + A^®n-92-.v) ^4' 

Equation (3) expresses that the added Mg""" ions substi¬ 
tuted for A-site Fe^"" and the Mn'*"' ions for B-site Fe^'^ 
without any change in the initial degree of inversion 
except j: = 0-6. The experimental integrated intensity ratios 
of the A and B sites, ^^Fe spectra and those calculated 
on the basis of cation distribution obtained from X-ray 
intensity data given in table 1, are in very good agreement 
with each other (table 2) for 0-2<x<0-9. 

The saturation magnetization (a) and the magnetization 
number (n^) at 300 K obtained from magnetization data 
for x = 0-1-0-9 are summarized in table 3. From field 
dependence of magnetization and observed magnetic 
moments at 300 K (table 3), it is clear that the samples 
with jc=0-0-0-6 show ferrimagnetic behaviour in agree¬ 
ment with Mossbauer data at 300 K (figure 3) which 
decreases with increasing x values and samples with 
x>0-7 display paramagnetic behaviour. Figure 4 shows 
the variation of with x for x = 0-0-0-9 at 300 K, 
increases with x up to x < 0-3 and thereafter it decreases. 

According to Neel’s (1950), two sublattice models of 
ferrimagnetism, the magnetic moment per formula unit 
in f. 1 ^, « 3 , is expressed as 

nlix) = M^{x)-M^{x). (4) 

where and are B and A sublattice magnetic 
moments in The nj (//.;^) values for x = 0-0-0-9 were 
calculated using the ionic magnetic moments of Fe^"", 
Mn'^'^ and Mg^'" with respective values S/ig, 3;rg and 


Table 3. Saturation magnetization (emu/g) and saturation 
magnetization per formula unit in Bohr magnetons, at 300 K. 
The values of the exchange constants in K used in calculating 
ctyK given as footnote. 



Magnetization (300 K) 

Yafet-Kittel angles 

X 

(emu/gm) 


«yk 

«YK 

0-0 

23-2 

0-80 

0° 

0° 

0-1 

29-2 

1-03 

0° 

0° 

0-2 

43-2 

1-50 

0° 

0° 

0-3 

45-3 

1-55 

14°36' 

19°40' 

0-4 

39-5 

1-32 

3r27' 

30°05' 

0-5 

38-4 

1-26 

40°09' 

39°22' 

0-6 

28-5 

0-92 

53°07' 

48°03' 

0-7 

22-2 

0-70 

63°ir 

57°44' 

0-8 

11-4 

0-36 

83°39' 

7r06' 

0-9 

- 

- 

- 

- 

4 = - 

14-0; J^ = - 

13-0; J, = - 

-10-0; J; = -4-5 

and =-7-5. 

OFb 

in (3) and 

the results 

are shown in 

figure 4 as a 

solid 

line. The 

calculated 

ng values for 

X = 0-0~0-3 at 


300 K are in good agreement with the experimentally 
found values confirming the collinear spin ordering, 
while for x = 0-4-0-9 they clearly differ from the observed 
values (figure 4) indicating that significant canting exists 
on B-sites, suggesting magnetic structure to be non- 
collinear. Thus the change of spin ordering from collinear 
to non-collinear displays a strong influence on the vari¬ 
ation of the saturation magnetic moment per molecule 




as observed by magnetization (figure 4) with chemical 
composition. 

In order to confirm the presence of a canted spin 
(non-collinear) structure indirectly, we have recorded 
Mossbauer spectra at 12 K for typical samples jc = 0-2, 
0-6, 0-8 and 0-9 as shown in figure 5 in addition to 
Mossbauer spectra taken at 300 K (figure 3). The spectra 
exhibit the coexistence of magnetic as well as very small 
paramagnetic component. The full lines in figure 5 
correspond to the least squares fitted spectra using two 
magnetic sextets and a central doublet. The magnetic 
components, which are well resolved at these tempera¬ 
tures, have been assigned to magnetically ordered Fe^"^ 
ions at the tetrahedral (A) and octahedral (B) sites. The 



V«locity(mm/s) 



hyperfine fields along with isomer shifts obtained from 
Mossbauer data analysis are listed in table 2 as a function 
of X for both A and B sites. Both hyperfine fields 
and //g) decrease with increasing x demonstrating reduc¬ 
tion in ferrimagnetic behaviour and magnetic coupling 
with increasing x values in agreement with mag¬ 
netization results (table 3). Isomer shifts IS (A) and IS (B) 
show very little change with Mn-content. The quadrupole 
shift is absent for all Mn-contents. 

An indirect check for the presence of a canted spin 
(non-collinear) structure follows from the apparent pro¬ 
portionality between II and the average sublattice 
magnetization. 

It can be shown that 
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I H^ {x)\ 

R 


M^{x)~ 


l\(x)l 

A 




(5) 


where /t(x) is the magnetic moment per formula unit 
(3), and I (x) I and I (;r) I are magnitudes of the 
average nuclear magnetic"' fields for the B and A sites, 
respectively. It is assumed that the relative saturation 
magnetization at 300 K is equal to ju(x)/ju(0) at 12 K. 
The values of magnetic moment per formula unit 
a function of a: were determined from 
(5) using the values of nuclear magnetic fields at 12 K 
from Mossbauer results (table 2) and the same are 
displayed in figure 4 and table 2. There is a reasonable 
agreement amongst the values of obtained from 

Mossbauer data, nj^(x) from magnetization measurements 
and ng(A:) calculated from Neel’s model as shown in 
figure 4 for a: = 0-0 to 0-2 confirming collinear magnetic 
structure. Mossbauer data, ng(x), and ^bCa:) for x>0-3 
clearly differ from magnetization data, WbCa:), indicating 



Figure 4. Vkiation of (Ug) with Mn concentration at 300 K 
for Mgj^^Mn^FCj. 2 ^, 04 . The dots with circles represent the 
magnetic moment obtained from Mossbauer data and the 
triangles represent data from magnetization measurements. The 
solid line corresponds to Neel’s model. 


that significant canting exists on B-site suggesting mag¬ 
netic structure to be non-collinear. 

Adopting the procedure of Kulkarni and Patil (1982), 
the condition for the existence of YK angle can be 
found for a: >0-3, and the molecular field acting on 
various ions in the present system are given by 


HA(Fe) 


2aa « a 

//g/Mn) 


a y' e' y s 


= 

p e' d' 8 d 

7/B2(Mn) 


ay s ye' 

/7g2(Fe) 


J3 £ d s' d' 


X 


(0-92 - a:) m^(Fe) 
l/2(A:)mm(Mn) 

1/2(108 -a:) mg,(Fe) 
l/2(x)mB2(Mn) 

1/2(1 08 - x) mg 2 (Fe) , 


( 6 ) 


where H^(Fe) is the molecular field acting on the Fe^'’’ 
ion on the A-site due to all other ions, etc and 


1 m^(Fe) 1 = I mg,(Fe) 1 = 1 mg 2 (Fe) 1 = 
I mB,(Mn) 1 = I mB2(Mn) I = S/^g, 
mg, • mA = -1 /Wg] M mA I cos 
rngj • ~ ^ I I '^A ^ ^YK ’ 

mg, • mgi = 1 Wg, 11 mB 2 1 cos . 


The molecular field constants are related to the fol¬ 
lowing sets of interactions 

= A(Fe) - A(F&), 
a = A(Fe) - Bj(Mn), 

^ = A(Fe)-B,(Fe), 
y = B,(Mn)-B"(Mn). 
y' = Bi(Mn)-B,(Mn), 

(5 = B,(Fe)-B 2 (Fe), 

(5' = B,(Fe)-B,(Fe), 
a=B,(Fe)~B2(Mn), 

£' = B,(Fe)-B,(Mn), 


where i = 1,2. 

The part of the interaction energy involving YK angles 
may be written as 

E(YK) = [ 15 ( 0-92 - x)xa + 25 ( 0-92 - jc) (1 -08 - jc) 

X cos 1 /4 [ 9 x^ y -f 25(1 -08 - a:)^ (5 + 30 a: 

X ( 1-08 - a) s] cos 2ayy^. 


(7) 



The energy is minimum for sin = 0, which corresponds 
to Neel’s configuration or for 

15X0-92 -x)a + 25(0-92 - x) (1 -08 - x)^ 
cn^ {y. = ——~~—■—--- 

9Xy +25(1-08-jc)-6 +30X1-08’ 

( 8 ) 

which corresponds to the situation where the YK ordering 
is feasible. Here, for example the molecular field constant 
is related to the exchange integral by 

2Z^3y[A(Fe)-B.(Mn)] 

^ \T 2 ’ 

<?A 



where the number of B; nearest neighbour to A, 
A^g, the number of ions per formula unit, p^, the Bohr 
magneton and g the Lande’s splitting factor. 

On physical grounds, it is expected that y', J^, J' 
will be close to 7^,, and J^, respectively, but they are 
not necessarily equal. 

The molecular field constants can be obtained from 
the observed variation of the saturation magnetization 
with Mn concentration. The exchange constants so 
obtained should be consistent with the existing data on 
similar ferrimagnetic systems. Kulkarni and Patil (1982) 
have found y = -5-25K for the CU|_ynye 204 system, 
in the case of Cu, _ynye 204 , is the interaction between 
A(Fe-^y and B.(Cu^'^), in the present case it corresponds 



Figure 5. Typical Mossbauer spectra of Mg, . Mn Fe,_-O, taken at 12 K for (a) x = 0-2, 0-6 
and (b) 0-8, 0-9. - - ^ 4 
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to the interaction between A(Fe^^) and Bj(Mn'’^). Assuming 
is approximately equal to the product of the moments 
of the ions, one gets - Cu^^) = - 5 K and y^CFe^"" 

-Mn''"') = -15 K. The value of 2^ = -5K compares rea¬ 
sonably well with J^ = -5-25 K value of Kuljcami and Patil 
(1982) for Cu, __,Zn/e 204 . So = - 14 K and 7^ = - 13 K 
values are accepted fpr Mg,^^,Mn^Fe 2 _ 2 ^P 4 . In a similar 
fashion the values of T,, J^, and have been estimated as 
7^ = - 10 K, J^, = - 4-5 K and 7, = - 7 5 K. The values of the 
exchange constants used in calculating Yafet-Kittel angles, 

from (8) are listed in table 3 along with values. 

The experimental values of have been obtained 
from the observed variation of n^ with x using the 
relation 

n= cos a^^-Mp). (9) 

These values have also been given in table 3. There is 
a satisfactory agreement between the experimental and 
theoretical values of Thus the observed variations 
of the saturation magnetization has been explained on 
the basis of the existence of Yafet-Kittel angles on the 
B-site spins. This suggests that A-B and B-B super¬ 
exchange interactions are comparable in strength. 
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Fracture toughness correlation with microstructure and 
other mechanical properties in near-eutectoid steel 
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Abstract. The variation of yield strength and fracture toughness was investigated for four different heat 
treatments attempted on specimens of a near-eutectoid steel. The aim of this study was to optimize the 
microstructure for simultaneous improvements in strength and toughness. Further, the fracture toughness 
deduced through empirical relations from tensile and charpy impact tests was compared with those measured 
directly according to ASTM Designation: E399. Among the four different heat treatments attempted in this 
study, the plane strain condition was valid in the fracture toughness tests for (i) normalized and (ii) hardened 
and tempered (500°C for 1 h) treatments only. The latter of the two heat treatments resulted in simultaneous 
improvement of strength and plane strain fracture toughness. The finely-dispersed carbides seem to arrest 
the crack propagation and also increase the strength. The pearlitic microstructure of the former leads to 
easy crack propagation along cementite platelets and/or cementite/ferrite interfaces. The nature of variation 
of empirically determined toughness values from tensile tests for different heat treatments is similar to that 
measured directly through fracture toughness tests, although the two sets of values do not match quantitatively. 
On the other hand, the toughness data deduced from charpy impact test is in close agreement with that 
evaluated directly from fracture toughness tests. 

Keywords. Fracture toughness; near-eutectoid steel; mechanical properties; heat treatment. 


1. Introduction 

One of the requirements of designing engineering struc¬ 
tures is to enhance the load-bearing capacity without 
plastic yielding of the components. At the outset, it may 
appear that an increase in yield stress would meet this 
requirement. But it is not so because of increasing danger 
of brittle fracture in high strength materials. Simultaneous 
improvements in yield strength and fracture toughness 
are needed rather than either of them. While continuum 
mechanics is the basis of structural design, the control 
of mechanical properties of structural materials is in the 
realm of materials science requiring an understanding of 
the structure-property relations. Fracture toughness being 
a more complex property than the yield strength, the 
crucial issue requiring understanding is the effect of 
microstructure on crack propagation. This will provide 
solutions to problems like: (i) a systematic development 
of tougher alloys (alloy design), (ii) the inverse rela¬ 
tionship between strength and toughness of a material. 


(iii) effect of heat treatment on toughness, and (iv) 
non-unique relation between strength and toughness. The 
variation in strength and fracture toughness for different 
microstructures in a near-eutectoid steel has been explored 
in this study. Further, the correlation between the fracture 
toughness measured directly and that deduced indirectly 
from other mechanical properties (obtained from tensile 
and impact tests) was also examined. The simpler test 
results may be availed to assess fracture toughness 
through such correlations. 

2. Experimental 

2.1 Materials and heat treatment 

The chemical composition of the near-eutectoid steel 
used in this study is shown in table 1. Starting with 
hot-rolled plate of 12-5 mm thickness, various specimens 
for charpy, tensile and fracture toughness tests were 
machined according to the dimensions specified in figures 


Table 1. Chemical composition of near-eutectoid steel. 

Element C Si Mn S P Fe 
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la-c. These specimens were coated with a ceramic slurry 
prior to heat treatment in order to minimize oxidation 
and decarburization. Four different heat treatments were 
attempted on these specimens in such a way that their 
strength and toughness varied with some overlap in their 
range. The heat treatment schedules for this purpose 
were selected from ASM handbook (from the listed 
properties for different heat treatments). The details of 
the adopted heat treatments are mentioned in table 2. 

2.2 Fracture toughness test 

The fracture toughness tests were performed according 
to ASTM standard E399-83. Fafigue precracking and 
fracture toughness tests were carried out on an MTS 
810 machine using three point bend test specimens. 

Fatigue precracking was carried out by cyclically load¬ 
ing the notched specimens at a load ratio of 0-1 in load 
control mode. The maximum load of the fatigue cycle 
was chosen in such a way that the maximum stress 
intensity factor did not exceed 80% of the anticipated 
plane strain fracture toughness {K^^) in each case. The 
precracked specimens were loaded monotonically in stroke 
controlled mode till fracture. From the recorded load 



( All dimensions, in mm ) 

Figure 1. (a) Charpy test specimen, (b) tensile test specimen 
and (c) three-point bend test specimen. 


versus crack opening displacement (COD) data, the 
apparent fracture toughness (ATq) was calculated using 
the relation: 

K^ = [P^S/BW^'^]f{a/W), ( 1 ) 

where is the load corresponding to 2% crack extension 
(obtained through the intersection of 5% secant line), S 
the span length, B the thickness, W the specimen width, 
a the crack length and/(a/W') is the specimen geometry 
and crack length related correction factor (known as Y 
factor), which is obtained analytically from the appropriate 
relations (ASTM E399; Provan 1988). The validity of 
plane strain condition is assessed by comparing the 
specimen thickness with 2-5 where is the 

0-2% yield stress in uniaxial tension. 

2.3 Charpy and tensile tests 

The impact energy (energy/unit area) of the standard 
charpy impact test is converted to equivalent making 

use of an empirical relation (Rolfe and Novak 1970): 

0-22 (CX\ (2) 

where E is the elastic modulus in MPa, the plane 
strain fracture toughness in the charpy energy 

in MPa.m and the units of the constant are [l/MPa.m]"’^. 
Standard tensile tests were also performed on differently 
heat treated specimens using an Instron testing machine. 
Equivalent values were deduced through an empirical 
relation (Hahn and Rosenfield 1968): 

K^^ = n(2Ea^^E/3y^\ (3) 

where n is the strain hardening exponent and the 

fracture strain in uniaxial tension. 

3. Results and discussion 

3,1 Strength and fracture toughness 

The tensile stress-strain curves and typical load- 
displacement (COD) curves of precracked specimens in 
fracture toughness tests are presented in figures 2 and 


Table 2. Tensile test data of heat treated specimens. 


SI 

no. 

Heat treatment"^ 

0-2% YS 
(MPa) 

UTS 

(MPa) 

Elong. 

(%) 

RA 

(%) 

Equivalent 
fracture toughness 
(MPa 

1. 

Normalized 

700 

1060 

19-1 

36-0 

30-2 

2. 

Oil hardened and tempered 
at 680°C for 1 h 

700 

780 

23-8 

58-4 

41-2 

3. 

Furnace cooled 

520 

912 

23-8 

36-7 

34-2 

4. 

Oil hardened and tempered 
at 500°C for 1 h 

1464 

1568 

7-2 

16-7 

38-7 


*Obtained from eauation f31: '^all austenitized at 860°C. 
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3, respectively. The tensile yield strength and fracture 
toughness (K^) data of differently heat treated specimens 
are given in tables 2 and 3. It is seen that the yield 
strength of normalized specimen and that tempered at 
680°C for I h is the same. The specimen tempered at 
500°C for 1 h exhibited the highest yield strength among 
all the heat treated specimens. It is well known that in 


a hardened steel, when heated successively to higher 
temperatures above room temperature, submicroscopic 
particles of cementite are formed through the gradual 
loss of excess carbon in martensite. With continued 
heating, the particles of cementite grow in size by a 
process of coalescence and growth. This coalescence is 
more pronounced at higher tempering temperature. The 


Table 3. Fracture toughness test data of heat treated specimens. 


SI 

no. 

Heat treatment^ 

Crack length 
a (mm) 

1. 

Normalized 

12-02 

2. 

Oil hardened and tempered 
at 680°C for 1 h 

12-08 

3. 

Furnace cooled 

11-96 

4. 

Oil hardened and tempered 
at 500°C for 1 h 

12-06 


’ IP 

max''' Q 

Fracture 
toughness (A'q) 
(MPaVm) 

0-2% YS (MPa) 
from tensile test 

1-10 

47-5 

700 

1-09 

84-7 

700 

1-00 

55-2 

520 

1-01 

69-0 

1464 















dispersion strengthening effect of particles depends on 
their size and spacing. The coarsening of cementite 
particles causes a softening of the steel at higher tempering 
temperature as seen in the case of treatment 2 (tempering 
at 680°C for 1 h). On the other hand, higher strength 
obtained for treatment 4 (tempering at 500°C for 1 h) 
results from the ultrafine cementite particles. 

The true stress-true strain curves of samples subjected 
to different heat treatments are shown in figure 2 and 
it may be noted that the true stress-true strain values 
of the last point of these stress-strain curves are based 



on the percent reduction in area at fracture. On checking 
the validity of plane strain condition in fracture toughness 
tests, the plane strain condition is valid only for the 
normalized (treatment 1) and tempered at 500°C for 1 h 
(treatment 4) cases. For fracture toughness tests involving 
the other two heat treatments, the condition of plane 
strain is not valid. Thus quantitative comparison of 
fracture toughness is possible for treatments 1 and 4 
only but not for 2 and 3. Considering the treatments 1 
and 4, the yield strength of the tempered (500°C) specimen 
is nearly twice that of the normalized one. It can be 
noted that despite the increase in strength, the fracture 
toughness of tempered (500°C for 1 h) specimen is 
significantly higher than that of normalized one. Thus 
the tempered (500°C for 1 h) specimen exhibits simul¬ 
taneous increase in strength and toughness relative to 
the normalized treatment. The usual inverse relation 
between strength and toughness is not applicable while 
comparing the properties arising from these heat 
treatments. The higher toughness of the tempered (500°C 



Figure 3. Typical load-displacement (COD) curves for near- Figure 4. Typical SEM fractographs for near-eutectoid steel: 
eutectoid steel: (a) Normalized and (b) hardened and tempered (a) Normalized and (b) hardened and tempered at 680°C for 
at 680°C for 1 h. 1 h. 
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Table 4. Charpy test data of heat treated specimens. 


SI 



Charpy energy (C^) 

Equivalent 
(MPa fm) 

no. 

Heat treatment 


(J/cm~) 

(MPa-m) 

1. 

Normalized 


13-6 

0-136 

48 

2. 

Oil hardened and 
at 680°C for I h 

tempered 

44-7 

0-447 

117 

3. 

Furnace cooled 


16-3 

0-163 

55 

4. 

Oil hardened and 
af 500°C for 1 h 

tempered 

21-7 

0-217 

68 


^Obtained from equation (2); '^all austenitized at 860°C. 


for 1 h) specimen follows from the tempered martensitic 
structure. The carbides are expected to be mostly spherical 
in shape and they arrest the crack propagation through 
their fine distribution. On the other hand, the normalized 
treatment results in pearlitic structure. In this case the 
cementite platelets and/or cementite-ferrite interfaces 
provide easy path for crack propagation (Curry and Knott 
1978; Porter et al 1978; Rama Rao 1978). Thus the 
tempered martensitic structure provides strengthening 
through the presence of ultrafine particles of cementite 
in ferrite and better toughness because the finely-dispersed 
carbides arrest crack propagation, 

Fractographic observations of normalized and furnace 
cooled specimens indicate the presence of quasi-cleavage 
facets and dimples (figure 4a). A change in plane of 
fracture from one grain to the other can also be noted 
in the cleavage dominated regions. The path of the main 
crack might have been the cementite platelets or the 
interface between cementite and ferrite lamellae of the 
pearlitic microstructure. These easy paths of cleavage 
fracture lead to low value of fracture toughness. 

In the tempered steels, fracture occurred predominantly 
by nucleation, growth and coalescence of voids. Micro¬ 
voids are usually initiated at oxide and sulphide inclu¬ 
sions, carbides and other imperfections. The carbide 
particles assist the void formation by providing weak 
interfaces which separate in tensile loading. The particles 
may also fracture and the fractured parts move apart by 
plastic flow of the matrix. The individual voids thus 
formed join .together by local necking of ligaments 
between the voids and there is some plastic flow before 
fracture occurs. This process is associated with more 
energy absorption and, thus, a higher toughness. 

3.2 Correlation between fracture toughness and 
other mechanical properties 

The estimates of from the tensile and impact test 
data are shown in tables 2 and 4, respectively. Considering 
the toughness estimates from tensile tests, one to one 


correlation is not seen between them and the apparent 
toughness {Kf) measured directly. However, the trends 
of change in toughness with heat treatment are similar. 
On comparing the data obtained from fracture tough¬ 
ness tests with those estimated from the impact energy 
values, good quantitative agreement can be noted between 
them. 

4. Conclusions 

(I) From a comparison between yield strength and fracture 
toughness of specimens with similar microstructures, the 
toughness is seen to decrease with the increase in strength. 

(II) The hardened and tempered (500°C for 1 h) treatment 
results in two-fold increase in strength with a significant 
increase in fracture toughness relative to the normalized 
one. The fine particles of cementite in ferrite cause 
strengthening and also arrest the crack propagation and 
thereby enhance the fracture toughness. On the other 
hand, the pearlitic microstructure is of lower strength 
and it provides easy paths for crack propagation along 
the cementite platelets and/or cementite/ferrite interfaces. 

(III) The variation in fracture toughness values deduced 
indirectly from tensile tests follow the same qualitative 
trend as those directly measured for different heat 
treatments, but the two sets of values differ from each 
other. On the other hand, the fracture toughness data 
obtained from charpy tests is in closer agreement with 
that measured directly for different heat treatments. 
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Abstract. Surface modification of metals and alloys using lasers has acquired considerable importance 
because the changes achieved are highly localized and superficial. In this process the mechanical and 
structural properties of the bulk are retained as such because of the high temperature gradient and high 
rate of change of temperature that are unattainable by conventional methods. The present report is on 
surface hardening of oil-hardened nonshrinking (OHNS) steel using Ruby laser in normal mode. It is found 
that surface hardening increases thrice the size of base material and the microstructure of the laser-heated 
zone changes from pearlite phase to martensite phase. 

Keywords. OHNS steel; ruby laser; surface hardening; material processing; Pearlite phase; Martensite 
phase. 


1. Introduction 

Heat treatment of materials by lasers gains popularity 
because specific region with pre-determined phase 
changes could be achieved. Various authors have studied 
surface hardening of steels of different compositions and 
grades, employing continuously working (CW) and pulsed 
lasers (Avner 1974; Gnanamuthu 1979; Hugel 1988; 
Lampman 1981). Markarov et al (1989) studied the 
influence of laser and electron beam treatments on the 
structure transformations and abrasive wear resistance of 
hardened steel 20KHN3A which is widely used in the 
manufacturing of drilling bits and other components 
.operating under conditions of high contact loads using 
CW COj laser. Bradley (1988) theoretically studied the 
hardened depth obtained by scanning steel surfaces with 
a laser beam. The calculated depth of hardening was 
found to be in agreement with experimental results for 
three steels processed with a CO, laser. Laser surface 
hardening of ferrous materials is a process used to 
improve the mechanical properties of highly stressed 
machine and tool parts. In this process only a thin layer 
of the surface is heated to austenization temperature 
prior to quenching leaving the interior of the work piece 
unaffected. If the power density of the laser beam is 
high, heat is generated at the surface at a rate higher 
than the heat conducted to the interior parts of the 
material and within a short duration the thin surface 


layer will reach austenizing temperature whereas the 
interior of the workpiece will be at its initial temperature. 
Laser transformation hardening of steel is currently em¬ 
ployed in many automobile and tool industries because 
of its technical advantages over the flame and induction 
hardening. It gives minimum distortion, ability to harden 
localized and inaccessible areas, high-speed processing 
and improved mechanical properties such as hardness, 
fatigue and wear resistance. Oil-hardened nonshrinking 
(OHNS) steel which comes under tool steel is widely 
used in manufacturing tools for machining and gauging 
such as machine screw taps, threading dies, chasers and 
milling cutters. Studying the surface hardening of this 
steel using ruby laser will be useful for manufacturing 
industries. The present report is on surface hardening of 
OHNS steel using a Ruby laser in normal mode. 

2. Experimental 

The material chosen for the study was OHNS steel, 
which covers under dimensionally stable steels. The 
percentage of carbon and its alloys in OHNS steel are: 
carbon, 0-95%; chromium, 0-5%; V, 0T5%; W, 0-6% 
and the rest is iron. The properties of OHNS steels are 
given in table 1. The microstructure of the specimen 
(parent material) is pearlitic and its microhardness 
280 VHN. This material is useful for threading tools of 
all types, such as machine screw taps, threading dies, 
milling cutters, wood working tools and gauges. 
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The specimen were cut into blocks of 5 x 3 x 1 cm 
and the roughly cut surfaces of the specimen are ground 
in the grinding wheel. The surface was polished by using 
polishing papers and alumina. The papers used for poli¬ 
shing are given in the increasing order. Rough polishing, 
120, 220, 400, 600; fine polishing, 1/0, 2/0, 3/0, 4/0, 
and alumina polishing is done on the specimen using 
rotary wheel. The sample is etched using ferric chloride 
to obtain the clear grain boundaries. 

The hardness of the base material is measured using 
Shimadzu micro hardness tester (load 200 g) and found 
to be 280 VHN. The sample was heat treated by using 
Ruby laser (Laser applications, Inc., Model No. 933R 
3L-2) of wavelength 694-3 nm, pulse width 300/^s, output 
energy 1-5 J, beam divergence 1 mrad and one pulse per 
min. The specimen is fixed in an optical upright which 
has vertical, horizontal as well as front and back move¬ 
ments. The laser beam was focussed by a glass lens of 
focal length 30 cm and the depth of focus was so 
adjusted that there was intense heating, but below melting 
or vapourization. If the laser beam is focussed it causes 
pitting of the sample. To avoid this the distance between 
the specimen and the lens are adjusted by trial and error 
method. The specimen is coated with carbon soot to 
avoid reflection of the incident radiation. The optimum 
range of the coating thickness is found to be 0-05 mm 
to 0-1 mm. 2 mW He-Ne laser is used for alignment. 
The sample is irradiated by Ruby laser with one pulse 
per min, and the hardness was measured. The diameter 
of the irradiated beam on the specimen is 0-5 mm. 
Experiments were repeated seven times. The microhard¬ 
ness for different laser shots are measured. The hardness 
at the edges (0-2 mm from the centre) and at the centres 
for different shots are measured and these values are 
given in table 2. The values show that the hardness is 

Table 1. Properties of OHNS steel. 

Thermal conductivity (K) 24-27 Wm~' k“' 

Thermal diffusivity (k : l^pC) 6-687 X 10’*’ m^/s 
Emissivity 0-8 

Microhardness 280 VHN 

Density {p) 8670 kg/m^ 

Specific heat capacity (C) 418-6 J/kg/k 

Microstructure present Pearlite structure 


Table 2. Microhardness at different distances. 


SI. No. 

Microhardness (VHN) 

Edge 

Centre 

Edge 

1 

560 

840 

665 

2 

569 

927 

680 

3 

532 

946 

610 

4 

520 

908 

670 

5 

543 

892 

603 


maximum at the centre and decreases at the edges. It 
may be due to the gaussian distribution of the beam 
profile. The maximum value of microhardness is found 
to be 946 VHN which shows that the hardness is about 
three times more than the hardness of the base material. 
Different laser shots are superimposed on the same spot 
and the microhardness was measured and the values are 
shown in table 3. The hardness is increased by three 
times for the first shot and after that the hardness slowly 
increases. 


Table 3. Microhardness for number of superimposed laser 
shots. 


No. of shots (superimposed) 

Microhardness (VHN) 

Base metal 

280 

1 

790 

2 

820 

3 

840 

4- 

850 

5 

857 




Figure 1. Microstructure of a. unhardened and b. hardened 
zones. 
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Figure 2. a-b. Microstructure of hardened-unhardened zones. 

The laser-treated surface is observed through optical 
microscope with magnification 200 x and it was found 
that the hardened surface has martensite phase whereas 
the microstructure of the parent material is pearlite phase 
which are shown in figures la and b respectively. The 
interface between the hardened and unhardened zones 
are shown in figures 2a and b. 

The surface temperature, T of the specimen is calculated 
using the formula (Rykalin et al 1988): 

T=2EPkt/AK 

where, E is the emissivity, P the power, k the thermal 


diffusivity, t the interaction time of laser, A the area 
and K the thermal conductivity. 

The hardened depth is found by using the expression 
(Rykalin et al 1988): 

Z=iAdt/str 

a is the beam diameter, and are the melting and 
hardening temperatures, respectively. 

The surface temperature and hardened depth are cal¬ 
culated and the values are 872°C and 40/^m, respectively. 

3. Conclusion 

The microhardness of the base material is 280 VHN. 
The hardness obtained by furnace hardening is 680 VHN. 
By using ruby laser the maximum hardness achieved is 
946 VHN. The laser heat treatment increased the hardness 
of the specimen to a larger extent. The comparison of 
microstructure of the laser-heated and unheated region 
showed a clear difference of pearlite and martensite 
matrices and the transition between the hardened and 
unhardened region was also observed. This study has 
showed the technical feasibility of laser surface hardening 
treatment to improve the wear resistance of parts of ring 
gear and engine valves without affecting the interior 
parts. 
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Abstract. An experimental setup is developed for the measurement of dynamic magnetoelectric effect (ME) 
in polycrystalUne materials, using a time varying DC magnetic field on which an AC magnetic field is 
superimposed. The experimental data on ME on BigPeTijOj^ and a solid solution of (90%)BiFe03-(10%)BaTi03 
are obtained using this setup. The linear and higher order ME coefficients can be evaluated from the output 
voltage. The temperature variation of ME data gives additional information corroborating with the data on 
magnetization. 
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1. Introduction 

The magnetoelectric effect (ME) in an insulator is the 
electrical polarization on applying DC or AC magnetic 
field (ME)j^ or the magnetic polarization on the application 
of DC or AC electric field (ME)^, It occurs due to the 
interaction between the magnetic and electric dipoles. 
The condition for the realization of ME output is that 
there should exist electrical and magnetic ordering in 
the material at a given temperature and also that the 
electrical resistivity is reasonably high (Rado and Folen 
1961, 1962; Freeman and Schmid 1975; Siratori et al 
1992). 

There are basically two kinds of materials which 
exhibit magnetoelectric output, namely single phase 
materials and composite materials. The single phase 
materials have an ordered lattice having two transitions 
corresponding to the ferroelectric to paraelectric, the 
other ferro/ferri/antiferro to paramagnetic state. While in 
the two-phase composites, one phase is piezomagnetic 
and the other piezoelectric. In the single phase material, 
ME output arises due to the interaction between the 
magnetic and ferroelectric sublattices (Suryanarayana 
1994) while in the composite, ME output, a product 
property, is essentially due to the interaction between 
two participating phases (Van Suchtelen 1992; Surya¬ 
narayana et al 1992). Further, it is to be noted that in 
neither of the individual phases ME signal is observed. 
The synthesis of the materials is itself a challenging 
task since the impurity phases and interphases that may 
be formed in the composite and other defects, if any. 
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may also affect the signal. In the case of single phase 
materials, most of the known ME compounds have 
complicated crystal structures and show output at low 
temperatures. However, the best known ME materials 
like CrjOj and boracites were studied in the single crystal 
form (Astrov 1960; Folen et al 1961; Ascher et al 1966). 
In most other cases, the data on polycrystalline materials 
can still yield useful information to understand the struc¬ 
ture-property correlations. 

An important aspect of magnetoelectric phenomenon 
is the necessity to pole the samples both electrically and 
magnetically. A proper poling methodology would ensure 
dipole reorientations so that the sub-lattices (in the case 
of single phase material) or different constituent phases 
(in the case of composites) would effectively participate 
to yield reasonable polarization, to be detected as a 
voltage or charge output. 

Both linear and higher order ME coefficients are 
possible, which have their origin in the expansion of 
density of stored free enthalpy g (E, H, T) written in 
terms of a power series. The order of ME coefficient 
also depends on the symmetry to which the crystal 
belongs (Rivera 1994). The magnitude of ME output is 
very small but still well measurable by any of the 
following methods. 

2. Measurement techniques 

There are three methods of measuring ME output, (i) 
In the static method, ME output (charge or voltage) is 
measured as a function of increasing magnetic field using 
an electrometer having high input impedance. While 
poling, there is a possibility of charges to get accumulated 
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at the grain boundaries which move towards the electrodes 
during the measurement. The discharge of the charges 
needs quite a lot of time, before actually the final output 
is stabilized (Hanumaiah et al 1994; Singh et al 1994; 
Suryanarayana 1994; James et al 1997). The accumulated 
charges in some materials may lead to erroneous con¬ 
clusions as we inferred earlier in the case of BijFeTijO,^ 
(Singh et al 1994), (ii) for more precise data, quasi-static 
method is employed, which involves a time varying DC 
magnetic field and a high impedance electrometer. The 
output is measured as a function of time. This technique 
was used by Rivera (1994) in single crystal boracites 
and Ye et al (1994) in CrjB^OijCl. But this method 
cannot be employed for a polycrystalline material as the 
charge build up at the grain boundaries still affects the 
ME output and (iii) dynamic method. 

2.1 Dynamic method 

In the dynamic method, the output ME voltage is mea¬ 
sured with a bias AC magnetic field while keeping the 
quasi-static measurement in act i.e. the measurement is 
carried out with a time varying DC magnetic field in 
the presence of an AC field. The bias AC magnetic 
field employed will not allow the charges to move 
towards the electrodes since a suitable signal with an 
appropriate frequency is used where the polarity of the 
signal changes with time. In what follows, the authors 
describe the details of dynamic method of measuring 
ME output that has been set up. The basic equations 
necessary to evaluate ME coefficients from the voltage 
are given below (Rivera 1994). 

When a DC magnetic field is applied, ME output 
voltage (V), in a polycrystalline material showing 


second-order enect can be given by me toilowmg 
equation, 

where a and /? are the coefficients of linear and quadratic 
components of ME. Suppose that if an AC field /Zj, is 
superimposed over the DC field //„, then the effective 
field is 

77 = //„ 4-/Z(, sin ip)t). 

A lock-in-amplifier tuned to the frequency 
o) {(0 = 271 f) would measure the output ME signal. Then 
we can write, 

Vc.{a+I^H,)h,=a^{H,)h,. 

When is zero, it is possible to find out the linear 
term a and hence the second-order term /? is evaluated 
in the presence of DC magnetic field. When H,, is not 
zero, the above equation shows that a varies with /?, 
giving rise to a pseudo-linear coefficient «*(//„). 

2.2 Experimental 

For the above measurement, a Helmholtz coil is mounted 
on the pole pieces of a DC electromagnet, whose magnetic 
field can be made time varying by interfacing it to a 
computer. The time varying DC magnetic field can be 
achieved by a programmable stepper motor attached to 
the potentiometer of the DC magnetic power supply. 
The block diagram is given in figure 1. 

A stepper motor with 0-7 kg-cm torque and having 



AC Magnetoelectric effect (Block Diagram) 

Figure 1. Block diagram of dynamic magneioelectric experimental setup. 
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1-2° angle per rotation is used for the present purpose. 
A closed loop control system was achieved by measuring 
the output of the Hall probe in the Keithley 181 Nano¬ 
voltmeter and calibrating the Hall voltage with the field. 
Using the GPIB, it is possible to control the input power 
supply of the magnet using the combination of stepper 
motor and potentiometer. The AC Helmholtz coils each 
of 130 turns were mounted on the DC pole pieces. The 
ratio of the magnetic induction and the current through 
the coils is independent of the frequency. The coils are 
excited with a low-amplitude high-frequency (1-006 kHz) 
signal, using the internal oscillator of the lock-in-amplifier 
(SR-530). Since the signal current was very small to 
drive the coils, a power amplifier was used to amplify 
the current. The AC magnetic field thus generated was 
calculated knowing the current carried by the coils using 
a Keithley 196 DMM. The whole setup was properly 
shielded to avoid stray pickups. 

Keeping the DC magnetic field zero, with increasing 
AC field, the ME signal was recorded which gives the 
linear coefficient (a). Then by properly selecting a fixed 
value of AC magnetic field, the DC field was swept 
using the stepper motor through the computer. The 
induced ME signal in the sample was measured with 
the lock-in-amplifier (SR530). With known values of 
a, /i,), and V, the value of second-order coefficient 
0?) can be calculated. Since, the total measuring time 
to complete the experiment is only few minutes, no time 
is available for the charges to accumulate in the sample, 
as in the case of pure DC (ME)j^ measurements. The 
signal output in voltage represents the true value. The 
same experimental setup, was used to take a temperature 
scan of the ME signal at fixed AC and DC bias fields 
by inserting the sample in a dewar in between the pole 
pieces. The data on (ME)^ at different temperatures 
would be useful to identify whether there are any magnetic 
anomalies in the material (Ikeda et al 1994). 

Using the above experimental setup, AC (ME) mea¬ 
surements were carried out on the following materials, 
which have been studied earlier for their synthesis, 
characterization and DC(ME)j, output in our laboratory 
(Singh et al 1994; Singh 1996; Mahesh Kumar et al 
1998); (i) Bi 5 FeTi 30,5 and (ii) (90%)BiFeO3-(10%) 
BaTiOj. 

In addition to the data at room temperature (RT), 
measurements were also made at 77 K, The temperature 
was measured using a Lakeshore silicon diode sensor 
DT-470. The samples were electroded with silver paint 
on both the sides and properly annealed. Two copper 
leads were attached on both the sides of the sample 
using epoxy resin. Prior to the measurement, the samples 
were poled electrically (15 kV/cm at 110°C) and 
magnetically (3 kOe at RT). 

The AC field (h^) employed was 18 Oe at 1 -006 kHz. 
A floating crystal was used to measure ME output and 


the differential inputs of the lock-in-amplifier were used 
to measure the output signal with coaxial cables. Along 
with the signal, the phase shift was also recorded. To 
minimize the noise due to cabling, a pre-amplifier SR-550 
was included in the circuit. 

As pointed out by Rivera (1994b), this experimental 
setup needs no calibration, contrary to (ME)^ method. 
Nevertheless, we have reexamined the data on 
BijFeTijOij, which was investigated by Singh (1996). 

3. Results 

3.1 Bi,FeTi,0,, 

Figure 2 shows ME output of Bi 5 FeTi 30 , 5 . At room 
temperature, value of the output voltage decreased with 
increasing DC magnetic field (with AC field super¬ 
imposed) up to 6 kOe. The shape of the curve is identical 
with that obtained by Singh (1996). It appears that the 
saturation fields are still higher. At 77 K, while the 
nature of the output was similar to that at RT, with 



Figure 2. Variation of magnetoelcclric output with DC bias 
magnetic field for Bi 5 FeTi 30,5 at room temperature and 77 K. 
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increasing field, the difference between the low field 
and high field value is small. On reversal of the applied 
magnetic field, both the curves do not retrace the same 
path. The remanance at RT is 16-9 mV/cm, while at 
77 K the value is 49-8 mV/cm. The calculated values 
of a at RT and 77 K are 0-1 mV/cm-Oe and 2-8 mV/cm- 



Figure 3. a. Variation of magnetoelectric output with DC bias 
magnetic field for (90%)BiFe03-(10%)BaTi03 at room tempe- 


Oe, respectively. The value of linear coefficient increases 
with decreasing temperature. At RT the second order 
coefficient (J5) is 1-37 x 10"® mV/cm-Oe^ and that at 77 K 
is 5-62X 10-‘^mV/cm-Oel 

3.2 (90%)BiFeO^-(10%)BaTiOj 

Ismailzade et al (1981) have reported that BiFe 03 -BaTi 03 
forms a series of solid solutions with different structural 
modifications. It is rhombohedral between 100-67% of 
BiFe 03 , 67% to 7% cubic and from 7% tetragonal. They 
also reported magnetoelectric measurements at room 
temperature for the composition 95 %BiFe 03 - 5 %BaTi 03 . 
From the polarization vs magnetic field plot, they ob¬ 
served that the polarization increased with increasing 
field and at about 7 kOe, the output suddenly decreased. 
Recently, we have synthesized this system and reported 
its dielectric and electrical properties (Mahesh Kumar et 
al 1998). 

Figure 3a shows ME output of (90%)BiFeO3-(10%) 
BaTi 03 at RT. The value of the output increases with 
increasing DC field, reaches a maximum at ~ 5-5 kOe 
and later decreases. This field is lower than what 
Ismailzade et al (1981) have observed for ( 95 %)BiFe 03 - 
( 5 %)BaTi 03 . Thus the mole% of BaTiOj in BiFe 03 has 
an effect on ME output and also on the magnetic field 
required to observe the peak in ME vs field measurement. 

Figure 3b shows ME output of the sample at 77 K. 
ME output in this case peaks at ~ 6 kOe. However, this 
peak is less pronounced than that observed at room 
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Figure 5. Variation of magnetoelectric output with temperature 
for (90%)BiFe03-(10%)BaTi03. 


temperature. The maximum value of ME output obtained 
is ~ l-6mV/cm. The value of ME output has increased 
with decreasing temperature. 

Figure 4 shows the magnetization vs temperature of 
the compound (90%)BiFe03-(10%)BaTi03. It is seen that 
the sample shows ferromagnetic nature, but as the tempe¬ 
rature approaches RT, the magnetization does not 
approach zero. It indicates that there may be yet another 
magnetic anomaly after the cross over from ferromagnetic 
state with a around 255 K. Whereas there is a slope 
change around 140 K, which may be attributed to the 
magnetic disorder in the system which appears at low 
temperatures. 

Figure 5 shows the magnetoelectric output vs tempera¬ 
ture for the sample measured at = 18 0e and = 
3 kOe in the temperature range 77-300 K. Anomalies 
were observed at ~ 149 K and at 238 K corroborating 
with the slope changes in the magnetization data (figure 
4), which indicates that the present experimental setup 
is useful not only to measure the ME output but also 
to corroborate small variations/fluctuations in the mag¬ 
netization data. 

4. Conclusions 

A versatile dynamic technique to measure ME output 
has been setup. The technique can be used to perform 
a few additional experiments to yield the following data: 
(i) magnetoelectric output at room temperature and at 
different temperatures, (ii) from the data one can evaluate 
the linear and quadratic ME coefficients (a and /?), 
remanance and saturation fields, (iii) if the DC field is 


swept, both to positive and negative polarities, a butterfly 
loop can be traced between the applied field and output 
voltage, (iv) the phase angle of the output signal can 
be detected, which can be exploited for the possible 
device applications, (v) temperature scan of ME output 
can indicate possible magnetic transitions, and (vi) fre¬ 
quency scan at a particular temperature can also be 
made. A flat frequency response would be ideal for the 
material to be used as a transducer (Bracke and Van 
Vliet 1981). 
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Density of micro-quantity liquids by the method of rise of drops in 
immiscible liquids 
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Abstract. The conventional methods such as specific gravity bottle method, pyknometer method, and the 
Westphal-balance method or the capillary tube method where accurate weighing is a problem cannot be 
employed when liquids are available in micro-quantities. The method described may be employed to determine 
the density of micro-quantity liquids (even up to 0-5 //I) and it is found to be simple and rapid. The method 
also allows the analyst to retrieve the sample for further analysis. 

Keywords. Density; micro-quantity liquids; rise of liquid drops; immiscible liquid column. 


1. Introduction 

Studies on the rise of liquid drops through an immiscible 
liquid column have been carried out by several authors 
(Klee and Treybal 1956; Saffman 1956; Bhattacharya 
and Venkateswarlu 1957; Houghton et al 1957; Harmathy 
1960; Goldsmith and Mason 1962; Griffith 1962; Davis 
and Acrivos 1966; Zieminski and Raymond 1968; Grace 
et al 1976; Clift et al 1978). Srinivasan et al (1996) 
have previously dealt with the problem of rise of liquid 
drops with six variables, F, D, u, cr, // and p alone (F, 
drag force; D, diameter of the liquid drop, u, terminal 
velocity acquired by the liquid drop; >j, viscosity of the 
liquid in the column, p, density of the liquid drop; a, 
density of the liquid in the column), to arrive at an 
expression for the drag force, F acting on the drop in 
rise through the method of dimensions. It has been 
suggested that the simple expression for the constant 
quantity, S occurring in the process of simplification of 
the drag force expression, F may be used to determine 
the density, p of the liquid drop. 

The expression for the density, p obtained (Srinivasan 
et al 1996) is 

p=[{2cj + X)-[X{X + Ao)f^]ll, (1) 

where 

X = S^j]/{r/uf, ( 2 ) 

S:=[{r/u?'Ho-p)V[tr-p^'% (3) 

r being radius of the drop. 


*Author for correspondence 


The physical property, density assumes a position of 
great importance in the identification of pure liquids. 
When liquids are available in bulk quantities, conventional 
methods like specific gravity bottle method, pyknometer 
method, Westphal-balance method etc may be employed 
to determine this physical constant. When liquids are 
available in small quantities and when neither the con¬ 
ventional methods nor the capillary tube method (where 
weighing is a problem) is adoptable, a new method has 
to be developed to determine the density of the 
micro-quantity liquids. 

In the present paper, therefore, the method and the 
expression for p suggested earlier (Srinivasan et al 1996) 
have been employed to determine the density, p of liquid 
drops and the results obtained have been presented. 

2. Experimental 

2.1 Equipment and liquids used 

A long graduated glass cylinder (of internal diam. 0-05 
m and height 1-5 m) with a small side tube attached to 
it at the bottom and sealed with a rubber septum (figure 
1); stop watch (Racer), accurate to 0 05 sec; Hamilton 
precision syringe, accurate to 0-01 ill; Ethylene glycol 
and soap oil (Ranbaxy purified grade); Hexane, Heptane, 
MIBK, Xylene, Cyclohexane, Benzene, Toluene, Iso-amyl 
acetate (Fisher purified grade); Turpentine (Chemlab 
purified grade); Water, Petrol, Kerosene, Naphtha and 
Diesel (distilled); Palm oil. Groundnut oil and Gingely 
oil (KKR refined grade); Coconut oil (VVD refined 
grade); Castor oil (Richardson refined grade); Sandal¬ 
wood oil (Vibhav purified grade), were used (table 2). 
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ihe liquids viz. ethylene glycol, water, chlorobenzene 
and bromobenzene, which were immiscible with the 
corresponding drop liquids (table 2) were selected as 
column liquids and used for filling the cylindrical column. 
Liquid drops of known volume were gently injected at 
the bottom of the liquid column using a graduated 
Hamilton Precision micro-syringe. When the drop rised 
freely and vertically without oscillation, the terminal 
velocity, u was determined by observing the time, t 
required by the liquid drop of radius, r to cover the 



Figure 1. The set up used for conducting rising drop experi¬ 
ments. 


of diameter, u then me drop radius was obtained from 
r = where r-D/2. All the experiments were 

conducted at room temperature (25 °C). 

3. Results and discussion 

The data presented in table 1 for ten liquid drop-liquid 
pair systems having three data points each shows that 
rill is approximately constant. The other seventeen liquid 
drop-liquid pair systems (table 2) for which the experi¬ 
mental results have been obtained satisfy the same. The 
density of the column liquids and the drop liquids, 
viscosity of the column liquids, interfacial tension between 
the liquid drop and liquid column given in table 2 were 
determined by the specific gravity bottle method, Ostwald 
Viscometer and the method of drops, respectively. The 
values of and 5 are given in 

table 3. Table 4 furnishes the observed density of the 
liquids (refer table 2) and the density of liquids determined 
from (1) and the estimated error in percentage. It may 
be seen from table 1 that for liquid drops of different 
radii of a given liquid pair system, r/u is a constant. 
The value of S calculated from (3) has been found to 
be approximately constant for all the liquid drop- 
liquid pair systems (table 3) and its mean value is 
0-313366 m"' s^ (table 3). This mean value of S is also 
found to be approximately equal and agree to the ex¬ 
perimental mean value of S (0-3176 m”'s') predicted 
previously for the thirteen liquid drop-liquid pair systems 
(Srinivasan et al 1996). The density values of the liquid 
drop estimated from (1) (column 4, table 4) using the 
experimental mean value of S (0-313366 m“'s^; table 3), 
r/u (table 2) and i] (table 2) show that they are comparable 
with the observed density value of the liquids determined 
by the specific gravity bottle method (column 3, table 
4). The estimated error in percentage (column 5, table 
4) is found to be less than 0-5% between the observed 
and the calculated density values (table 4). 

Since by knowing the radius, r of the liquid drop and 
the terminal velocity, u attained by it in the liquid 
column, one may determine the density of the liquid 
drop from (1) with the experimentally predicted mean 
value of S (0-313366 m“‘s", table 3) and the viscosity 
of the liquid in the column (table 2), the only unknown 
quantity to be determined being r/ii. This may be 
accomplished in one or two min. 

The minimum amount of liquid sample required for 
this method is less than 1 /d, and therefore this method 
may be adopted to determine the density of a liquid 
sample available in small or micro-quantities for which 
the density cannot be determined by any other conven¬ 
tional method. It also provides a solution to the analyst 
who prefers, as far as possible to preserve the original 





Table 1, Experimental data for liquid drop-liquid pair systems. 
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Table 2. Liquid drop-column liquid pair systems and their physical constants. 


SI. 

no. 

Liquid 

drop 

Liquid column 

P. 

(kg-^) 

a 

(kg-^) 

(a-p) 

(kg-^) 

(Nsm ■) 

y , 

(nm ’) 

(X 10"^ ) 

r/u 

(s) 

1. 

Hexane 

EG* 

665-12 

1108-00 

442-88 

0-01520 

21-6 

0-017192 

2. 

Petrol 

EG* 

715-38 

1108-00 

392-62 

0-01520 

16-7 

0-018964 

3. 

Heptane 

EG* 

720-24 

1108-00 

387-76 

0-01520 

20-8 

0-019244 

4. 

Naphtha 

EG* 

733-29 

1108.00 

374-71 

0-01520 

15-3 

0-019708 

5. 

MIBK 

EG* 

792-85 

1108-00 

315-15 

0-01520 

3-1 

0-022752 

6. 

Kerosene 

EG* 

797-22 

1108-00 

310-78 

0-01520 

11-9 

0-022998 

7. 

Diesel 

EG* 

8.30-65 

1108-00 • 

277-35 

0-01520 

15-6 

0-025111 

8. 

Soap oil 

EG* 

857-01 

1108-00 

250-99 

0-01520 

18-0 

0-027218 

9. 

Xylene 

EG* 

857-50 

1108-00 

250-50 

0-01520 

13-9 

0-027318 

10. 

Benzene 

EG* 

870-60 

1108-00 

237-40 

0-01520 

9-1 

0-028376 

11. 

Palm oil 

EG* 

876-66 

1108-00 

231-34 

0-01520 

9-6 

0-028816 

12. 

Groundnut oil 

EG* 

910-91 

1108-00 

197-09 

0-01520 

12-0 

0-032461 

13. 

Gingely oil 

EG* 

915-85 

1108-00 

192-15 

0-01520 

12-3 

0-033150 

14. 

Coconut oil 

EG* 

917-27 

1108-00 

190-73 

0-01520 

13-3 

0-033273 

15. 

Castor oil 

EG* 

925-72 

1108-00 

182-28 

0-01520 

7.2 

0-034431 

16. 

Sandalwood oil 

EG* 

960-20 

1108-00 

147-80 

0-01520 

14-1 

0-040103 

17. 

Heptane 

Water 

720-37 

1000-00 

279-63 

0-00100 

38-7 

0-009765 

18. 

Cyclohexane 

Water 

775-04 

1000-00 

224-96 

0-00100 

20-4 

0-011578 

19. 

Kerosene 

Water 

797-34 

1000-00 

202-66 

0-00100 

43.8 

0-012480 

20. 

Soap oil 

Water 

857-15 

1000-00 

142-85 

0-00100 

34-2 

0-015956 

21. 

Xylene 

Water 

857-95 

1000-00 

142-05 

0-00100 

29-0 

0-016010 

22. 

Turpentine 

Water 

860-03 

1000-00 

139-97 

0-00100 

41-5 

0-016448 

23. 

Toluene 

Water 

860-89 

1000-00 

139-11 

0-00100 

39-7 

-0-016543 

24. 

Benzene 

Water 

870-78 

1000-00 

129-22 

0-00100 

35-0 

0-017382 

25. 

Iso-amylacetate 

Water 

882-15 

1000-00 

117-85 

0-00100 

29-2 

0-018537 

26. 

Water 

Chlorobenzene 

1000-00 

1097-99 

97-99 

0-00071 

46-1 

0-019610 

27. 

Ethylene glycol 

Bromobenzene 

1108-00 

1492-21 

384-21 

0-00085 

11-2 

0-008580 

P. ' 

Density of the liquid drop; a, density of the column liquid; 

j/, viscosity of the column liquid; y, interfacial 

tension between the liquid drop 


and the liquid in the column; /•, radius of the liquid drop; u, terminal velocity of the drop. 

(r/w): Mean experimental values. *Ethylene glycol. 

Table 3. The values of and p^^V(a-p) and S. 


SI. 


Liquid 


p^''V(a-p) 


no. 

Liquid drop 

column 

S 

1. 

Hexane 

EG* 

0-018284 

0-058232 

0-313985 

2. 

Petrol 

EG* 

0-021182 

0-068123 

0-310938 

3. 

Heptane 

EG* 

0-021653 

0-069211 

0-312855 

4. 

Naphtha 

EG* 

0-022441 

0-072267 

0-310529 

5. 

MIBK 

EG* 

0-027836 

0-089346 

0-311553 

6. 

Kerosene 

EG* 

0-028289 

0-090852 

0-311375 

7. 

Diesel 

EG* 

0-032276 

0-103915 

0-310600 

8. 

Soap oil 

EG* 

0-036422 

0-116637 

0-312268 

9. 

Xylene 

EG* 

0-036623 

0-116898 

0-313290 

10. 

Benzene 

EG* 

0-038771 

0-124287 

0-311947 

11. 

Palm oil 

EG* 

0-039676 

0-127986 

0-310003 

12. 

Groundnut oil 

EG* 

0-041431 

0-153134 

0-309774 

13. 

Gingely oil 

EG* 

0-048956 

0-157492 

0-310848 

14. 

Coconut oil 

EG* 

0-049228 

0-158792 

0-310016 

15. 

Castor oil 

EG* 

0-051821 

0-166917 

0-310460 

16. 

Sandalwood oil 

EG* 

0-065139 

0-209655 

0-310696 

17. 

Heptane 

Water 

0-030515 

0-095982 

0-317924 

18. 

Cyclohexane 

Water 

0-039396 

0-123753 

0-318344 

19. 

Kerosene 

Water 

0-044088 

0-139332 

0-316424 

20. 

Soap oil 

Water 

0-063736 

0-204950 

0-310983 

21. 

Xylene 

Water 

0-064060 

0-206200 

0-310669 

22. 

Turpentine 

Water 

0-066707 

0-209518 

0-318383 

23. 

Toluene 

Water 

p-067285 

0-210918 

0-319010 

24. 

Benzene 

Water 

0-072469 

0-228362 

0-317343 

25. 

Iso-amylacetate 

Water 

0-079810 

0-252024 

0-316676 

26. 

Water 

Chlorobenzene 

0-103059 

0-322714 

0-319351 

27. 

Ethylene glycol 

Bromobenzene 

0-027260 

0-086636 

Mean 

0-314650 

0-313366 


*Ethylene glycol; {r/u\ ?},p and (a-p) from table 2. 

/..^3/2 ,..t/2i yr _l/2 
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Table 4. Comparison of density values and estimated error in 
percentage. 


Density p (kg“^) 

-Estimated 

Estimated error 


no. 

Liquid drop 

Observed 

from (1) 

(%) 

I. 

Hexane 

665-12 

665-77 

- 0-098 

2. 

Petrol 

715-38 

712-98 

0-335 

3. 

Heptane 

720-24 

719-74 

0-069 

4. 

Naphtha 

733-29 

730-61 

0-365 

5. 

MIBK 

792-85 

791-32 

0-193 

6. 

Kerosene 

797-22 

795-56 

0-208 

7. 

Diesel 

830-65 

828-53 

0-255 

8. 

Soap oil 

857-01 

856-24 

0-090 

9. 

Xylene 

857-50 

857-44 

0-007 

10. 

Benzene 

870-60 

869-65 

0-109 

11. 

Palm oil 

876-66 

874-45 

0-254 

12. 

Groundnut oil 

910-91 

908-85 

0-226 

13. 

Gingely oil 

915-85 

914-44 

0-154 

14. 

Coconut oil 

917-27 

915-41 

0-203 

15. 

Castor oil 

925-72 

924-17 

0-167 

16. 

Sandalwood oil 

960-20 

959-02 

0-123 

17. 

Heptane 

720-37 

723-73 

-0-466 

18. 

Cyclohexane 

775-04 

778-12 

-0-397 

19. 

Kerosene 

797-34 

799-08 

-0-218 

20. 

Soap oil 

857-15 

856-14 

0-118 

21. 

Xylene 

857-95 

856-14 

0-133 

22. 

Turpentine 

860-03 

862-07 

-0-237 

23. 

Toluene 

860-89 

863-17 

-0-265 

24. 

Benzene 

870-78 

872-29 

-0-173 

25. 

Iso-amylacetate 

882-15 

883-31 

-0-131 

26. 

Water 

1000-00 

1001-75 

-0-175 

27. 

Ethylene glycol 

1108-00 

1109-33 

-0-120 


*From table 2. 

liquid samples in small quantities for identification and 
confirmation through other analytical means apart from 
density determination. 

Since the measurement of weight of the liquid drop 
is not involved in this method, the availability of a high 
precision balance is not an essential requirement as in 
the case of other conventional methods viz. specific 
gravity bottle method, pyknometer method, Westphal- 
balance method (Serence et al 1970; O’Hara and Oster- 
burg 1974; Giancoli 1984) etc for getting the accurate 
density value, and hence this method may be employed 
in any operational laboratory where high precision 
and expensive equipment is not available for accurate 
weighing. 


The added advantage of this method is that the liquid 
drop injected at the bottom can be retrieved from the 
top of the liquid column using a micro-syringe or filter 
paper for further analysis as the drop liquid and the 
column liquid are immiscible. 

4. Conclusion 

This method may be employed to determine the density 
of liquids available in small or micro-quantities and also 
provides a practical alternative to other conventional 
methods. The authors do not suggest that this approach 
is an acceptable substitute for those liquids which do 
not have any suitable column liquid. Therefore, the 
column liquid should be selected in such a way that it 
is immiscible with as well as denser than the liquid to 
be tested. 
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